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SECTION  1 


INTRODUCTION 


1. 1 Purpose  of  the  Study 

The  broad  objective  of  the  work  described  here  has  been  to  conduct  a 
design  study  aimed  at  developing  and  evaluating  techniques  for  suppressing 
false  alarms  in  air-defense  radars.  More  specifically,  the  output  of  the  study 
should  define,  to  the  extent  allowed  by  the  level  of  effort  expended,  circuitry 
required  in  a False  Target  Suppression  (FTS)  signal  processor  that  would 
implement  these  techniques.  This  processor  would  be  compatible  with  the 
U.S.  Army  Experimental  Array  Radar  (EAR)  which  is  presently  installed  at 
the  U.S.  Army  Missile  Research,  Development  and  Engineering  Laboratory 
at  Redstone  Arsenal,  Alabama. 

1.2  Organization  of  the  Study 

Pursuant  to  the  above  objective,  the  work  (as  well  as  this  report)  has 
been  organized  into  the  following  categories: 

i 

1.2.1  Background 

The  characteristics  of  the  EAR  and  of  its  existing  MTI  processor  have 
been  reviewed.  The  characteristics  of  the  target  and  clutter  environment 
within  which  it  can  be  expected  to  work  have  been  postulated.  The  charac- 
teristics of  other  processors  which  have  similar  FTS  functions,  both 
developed  within  Raytheon  or  available  in  the  literature,  have  been  reviewed. 

1.2.2  General  Processor  Configuration 

Based  on  the  background  information,  a general  FTS  block  diagram 
has  been  generated.  This  forms  a baseline  for  analytical  performance  studies 
as  well  as  detailed  processor  architecture  studies  and  the  cost/performance 
tradeoffs  which  follow. 


1.2.3  Analysis 

A detailed  analysis  has  been  carried  out  of  the  performance  of  the  FTS 
processor  (basically  expressed  as  target  detection  probability  for  a fixed 
false  alarm  rate)  as  a function  of  certain  parameters  such  as  number  of  target 
samples  coherently  processed,  input  signal-to-clutter  ratio,  input  signal-to- 
noise  ratio,  presence  or  absence  of  MTI  canceller  processing,  etc.  This 
study,  although  by  no  means  exhaustive,  does  clearly  indicate  the  performance 
improvements  that  might  be  anticipated  by  the  use  of  a properly  configured 
FTS  processor. 

1.2.4  Hardware/Software  Configurations  and  Tradeoffs 


Carrying  out  in  more  detail  the  work  outlined  in  1.  2.  2 above,  this 
task  examines  various  alternative  approaches  to  the  design  of  the  processor. 
More  detailed  block  diagrams  are  generated,  and  estimates  are  made  of 
such  cost -influencing  factors  as  the  number  of  components/cards  required  by 
each  approach  and  the  software  requirements,  as  a function  of  system 
specifications  such  as  number  of  instrumented  range  cells  and  number  of 
samples  coherently  processed. 

The  information  generated  in  1.2.4  is  combined  with  rough  estimates  of 
relative  development  costs  to  form  curves  of  the  cost  of  various  design 
approaches  as  a function  of  required  performance,  or,  conversely,  the 
performance  that  could  be  expected  for  a given  fixed  cost  investment.  This 
latter  is  perhaps  the  first  step  toward  the  application  of  "design-to-cost" 
philosophies  to  signal  processor  design. 


SECTION  2.  BACKGROUND  AND  ASSUMPTIONS 


The  Experimental  Array  Radar  (EAR 


The  EAR  ia  a C -band  phased  array  radar  designed  to  fulfill  both 
search  and  track  functions  in  air  defense  applications  (References  1 and  2). 
Two  crossed  beams- -vertical  and  horizontal  fans--are  generated  and  steered 
from  two  linear  arrays.  False  alarm  suppression  is  of  primary  significance 
only  in  the  search  function,  so  that  this  is  the  basic  function  considered  in 
this  study,  although  certain  of  the  processing  elements  described  here  such 
as  the  doppler  processor  are  also  of  considerable  value  in  improving 
tracking  performance. 

For  search,  only  the  vertical  fan  beam,  which  scans  in  azimuth  over 
approximately  £ 50°,  is  used.  Those  radar  characteristics  which  are 
important  to  the  signal  processor  are  listed: 

a.  0.2  microsecond  effective  pulse  width. 

b.  Nominal  PRF  of  5 kHz.  This  is  however  variable,  and  in  fact,  is 
almost  sure  to  be  varied  in  a programmed  fashion  to  overcome  target  blind 
speed  problems  and/or  as  an  ECCM  measure.  We  have  assumed  5 kHz  as 
an  upper  limit  on  the  PRF.  This  implies  a maximum  of  1000  range  cells  per 
PRF  interval  and  a 30  km  maximum  range,  although  in  practice  a smaller 
number  is  likely  to  be  instrumented --about  700--in  order  to  allow  some 
radar  dead  time  for  uncompressed  pulse  length  (if  pulse  compression  is  used), 
beam  steering,  etc. 

Q 

c.  2 azimuth  beamwidth,  implying  50  azimuth  angle  resolution  cells 
per  scan.  The  nominal  scan  time  is  assumed  to  be  about  0.  5 second  (or  10 
milliseconds  per  angle  dwell),  but  it  should  be  stressed  that  the  EAR  is  an 
experimental  radar  so  that  the  ability  to  vary  such  parameters  is  an  important 
part  of  the  design  of  any  processor. 


[ 


d.  There  is  a 3-pulse  digital  MT1  presently  implemented  with  the 
EAR.  This  is  a non-moving-window  design  with  one  output  per  three  pulses 
input.  The  FTS  processor  should  be  designed  so  that  some  appropriate 
elements  of  it  can  interface  with  this  MTI. 

2.  2 Limits  on  the  Processor  Design 

As  mentioned  above,  variability  is  an  important  property  of  the  FTS 
processor  to  be  designed.  But  to  keep  the  design  within  some  limits  we  must 
consider  certain  inherent  limits  on  the  radar's  output,  which  are  functions  of 
an  assumed  target  environment. 

Assume  that  targets  have  a maximum  velocity  of  Mach  3 (about  1000  ny^ec). 
The  minimum  velocity  is  to  be  as  low  as  possible.  The  clutter  environment, 
discussed  in  a later  section  is  typically,  fields,  wooded  hills,  weather,  etc. 

The  antenna  scan  itself  poses  no  limits  on  dwell  time  (and  therefore 
pulses  to  be  processed  within  one  beamwidth).  Considered  here  are  other 
limitations,  which  apply  largely  to  high-speed  targets  which  in  a tactical 
radar  are  likely  to  represent  the  greatest  threat.  These  are: 

• Scan  time:  The  rate  at  which  new  information  is  needed. 

• Range  walk:  The  motion  of  the  target  through  a range  cell  during 
the  time  on  target. 

• Target  acceleration  effects:  A limit  on  coherent  processing. 

• Maximum  target  coherence  time. 

2.2.1  Scan  Time  Limitation 

The  nominal  scan  time  is  0.  5 second,  corresponding  to  a dwell  time  of 

Q O 

10  msec  (for  2 beamwidth,  100  search  sector).  A more  generous  limit 
could  be  established  by  postulating  that  a Mach  3 target  be  detected  after 
penetrating  not  more  than  5 km  into  the  radar  coverage:  this  leads  to  a 
scan  time  of  5 seconds,  and  a dwell  time  of  100  msec.  This  latter  is 
probably  an  upper  bound  since  it  allows  no  time  for  radar  track  functions. 
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2.2.2  Range  Walk 

Given  a target  radial  velocity  o£  m/sec.,  a target  will  pass  through 
the  0.2  jusecond  (30  m)  range  cell  width  in  T^^  = 30/V^  seconds.  Beyond 
this  limit  there  is  no  integration  gain  within  a range  cell,  although  there  still 
will  be  an  increase  in  overall  target  detectability  due  to  detection  opportunity 
on  adjacent  range  cells.  For  a Mach  3 target,  this  limit  is  30  milliseconds 
(150  pulses  at  the  nominal  5 kHz  PRF). 


2.2.3  Target  Acceleration 

Here  we  assume  that  the  target  is  moving  in  a straight  line  at  constant 
velocity:  the  range  acceleration  is  due  to  the  rate  of  velocity  change  for 

such  targets  on  near -tangential  courses.  It  is  easy  to  show  (see  the  sketch) 
that  the  radial  acceleration 


At 


is  given  by: 


At 


.,2  2 . 

V cos  0 

R 


i 


Where  V is  the  target  velocity,  Rq  is  the  range,  0 is  the  off -tangent 
target  angle.  At  is  the  dwell  time,  and  AVR  is  the  change  in  radial  velocity 
during  the  dwell  time. 


Coherent  processing  will  not  be  efficient  if  AV^  equals  or  exceeds  the 
width  of  a "doppler  bin"  (measured  in  the  same  units)  of  the  doppler 
processor.  The  doppler  bin  width  fg  for  full  coherent  processing  during 
the  time  At  is  l/At  Hz.  The  equivalent  velocity  bin  width  is: 


A 


Then,  for  0 near  zero: 


< V O' 


2-3 


I 


f 


r 


i 


# 


i 

I » 

'+ 

l 

, f- 

> 


f 


For  A.  = 5.  5 x 10"2  m,  V = lO^m/sec.,  and  RQ  (minimum  range)  = 103m, 
At  should  not  exceed  5.2  msec.,  or  the  number  of  pulses  coherently  inte- 
grated (N)  need  not  exceed  26.  At  maximum  30  km  range.  At  is  28.  5 msec, 

N = 142. 

2.2.4  Target  Coherence 

From  Reference  3,  page  397,  the  doppler  spreading  due  to  target 

rotation  is  given  by: 

2 u > L 

l = -r- 

MAX  \ 

Where  uir  is  the  target's  real  or  apparent  rotation  rate  and  L is  the 
target  length.  For  a crossing  target  as  before,  the  apparent  target  rotation 
rate  is  given  by: 

v 

0)  = COS  0 

r Ro 

7T  V 

Then  fMAx  " fCT~  f° r a near  tangential  target, 
o 

3 3 

For  V = 10  m/sec,  a target  length  of  5 meters,  and  Rq  = 10  meters 

as  before, 

fMAX  = 180  Hz- 

The  corresponding  useful  period  of  target  coherence  is  At  = - — - — =5.5  msec, 

£MAX 

or  27  pulses.  At  maximum  range.  At  = 165  msec.,  N = 825. 

The  simple  relationships  derived  in  sections  2.2.3  and  2.2.4  are 
plotted  in  Figure  2-1.  These  should  not  be  interpreted  as  meaning  that 
longer  dwell  times  than  those  indicated  cannot  be  used,  but  only  that  the  full 
signal-to-noise  improvement  that  should  be  expected  from  coherent  inte- 
gration will  not  be  achieved  beyond  this  dwell  time.  Since  S/N  may  not  be 
of  importance  at  short  range,  the  dwell  time  limits  at  maximum  range  are 
probably  the  most  important  ones. 
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2.3  The  Clutter  Environment 
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2.  3.  1 ('» round  Clutter 

At  the  maximum  range  of  30  km  the  illuminated  ground  clutter  patch  is 
3 2 

30m  x 1050m  - 31  x 10  m in  area.  Using  C-band  data  from  Reference  4, 

page  262.  the  84th  percentile  clutter  coefficient  O is  >27  dB  for  wooded  hills.  ] 

3 ° 2 ] 

This  leads  to  a clutter  cross-section  of  31  x 10  x .002  = 62m  . Assuming 

2 

the  smallest  target  size  to  be  0.  lm  , the  distributed  C/S  ratio  is  about  30  dB.  1 


This  neglects  the  effect  of  point  clutter,  which  can  be  considerable. 
Cross-sections  of  10  m for  point  clutter  are  not  unusual  (Reference  5). 

This  implies  a C/S  ratio  of  50  dB  at  these  points.  Since  there  are  compara- 
tively few  of  such  points,  they  are  a logical  candidate  for  removal  by  clutter 
mapping. 

The  spectrum  of  the  distributed  ground  clutter  (from  Reference  4, 
page  274)  is  about  0.  5 m/sec.  in  width  (+  O)  under  very  strong  (~50  kt)  wind 
conditions,  proportionately  less  for  lighter  winds.  At  C-band  this  corres- 
ponds to  a doppler  band  of  18  Hz.  This  is  actually  + 9 Hz,  centered  at  zero 
frequency. 

2.  1.  2 Precipitation  Clutter 


The  volume  of  rain  illuminated  is  equal  to  the  ground  patch  area  of 

31  x I0^m  times  the  height  of  the  rainstorm  or  the  antenna  vertically- 

illuminated  altitude,  whichever  is  least.  For  an  assumed  storm  height  of 

8 3 

3 kin,  the  first  condition  governs,  and  the  volume  is  about  10  m at  maximum 
range. 


The  backscatter  cross-section  of  the  rain  is  the  reflectivity  coefficient 

at  C-band  times  the  volume.  The  reflectivity  coefficient  u,  from  Reference  3, 

page  106,  is  tabulated  below  for  various  rainfall  rates: 

o n.  a 2/3  C.at  max.  range 

Rain  Rate,  mm/hr.  u,  m /m  rain 


Light 

1 

5 x 10 

Moderate 

4 

7 x 10 

Moderate 

10 

2 x 10 

Heavy 

20 

io-6 

0.  5m ^ 
7.0 
20 
100 
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Typically,  (Reference  4,  page  211)  the  rain  clutter  has  a spectral 
width  of  up  to  5 m/eec.  about  a canter  velocity  of  + 30  m/sec. , depending 
of  course  on  wind  velocity  and  antenna  direction.  At  C~band  the  doppler 
width  is  about  180  Hz. 
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SECTION  3.  THE  GENERAL  APPROACH 


The  approach  described  here  is  a fully  digital  processor,  both  for 
compatibility  with  the  existing  EAR  equipment  and  because  analog  approaches 
are  not  nearly  powerful  or  flexible  enough.  Several  techniques  are  considered 
in  combination,  each  reinforcing  the  other  to  attack  the  variety  of  false 
target  sources  which  are  expected:  distributed  ground  clutter,  strong  point 
ground  clutter,  precipitation  clutter,  countermeasures  and  pulsed  interference. 

The  techniques  are: 

• Doppler  processing,  as  implemented  by  the  Fast  Fourier 
Transform  (FFT)  technique. 

• Range -cell  averaging  CFAR 

• Clutter  mapping 

• Frequency-cell  CFAR 


Although  it  is  recommended  that  all  of  these  techniques  be  implemented 
to  obtain  maximum  false  target  suppression  with  minimum  target  loss,  it 
is  not  necessary  that  they  all  be  implemented  in  order  that  benefit  be 
obtained.  For  example,  the  range-cell  CEAR  alone  when  used  with  the 
presently  existing  EAR  will  be  reasonably  effective. 

Figure  3-1  illustrates  the  proposed  configuration.  The  A/D  and  the 
lower  MTI  already  exist.  The  other  blocks  represent  new  equipment. 
Describing  in  general  terms  the  functions  of  each  block: 


The  existing  MTI  processes  input  data  in  block  form  producing  one 
output  pulse  for  each  three  input  pulses.  It  is  not  suitable  as  an  input  to 
the  doppler  processor.  However  it  can  be  input  to  a cell-averaging  CFAR 
processor,  and  it  is  shown  as  being  capable  of  being  switched  in  at  that  point. 


{ 
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General  Block  Diagram  of  EAR  Signal  Processor 


f 
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The  sliding  window  MTI  is  suitable  for  input  to  the  FFT,  and  is  so 
shown.  Although  "optional"  in  the  sense  that  the  FFT  itself  by  virtue  of  its 
doppler  filtering  behavior  has  certain  MTI  properties,  it  is  shown  in  Section  4 
that  the  implementation  of  an  MTI  at  this  point  will  enhance  the  total  per- 
formance of  the  FTS  processor  in  a clutter  environment. 

3.  2 Input  Weighting 

The  FFT  is  a block  processor,  processing  sequential  pulse  trains 
which  are,  most  efficiently,  powers  of  two  in  length.  The  truncation  implicit 
in  the  processing  of  a finite  block  like  this  will  give  rise  to  sin  x/x  side  lobes 
in  the  frequency  domain,  with  the  worst  sidelobe  some  13  dB  down.  By  pre- 
multiplying the  input  pulse  sequence  by  some  sort  of  a smooth  amplitude 
weighting  function,  the  sidelobe  level  can  be  controlled  to  any  desired  degree, 
with  a price  being  paid  in  widening  the  bandwidth  of  the  doppler  filter  and  in 
reduced  signal-to-noise  ratio.  The  Taylor  weighting  function  is  the  most 
efficient  and  is  the  one  considered  here. 

3,  3 Input  Buffer 

This  is  actually  a double  buffer.  A complete  block  of  input  data 
(N  pulses  by  M range  cells)  is  fed  in  real  time  into  one  buffer  while  the  FFT 
ia  processing  the  contents  of  the  other  buffer,  which  had  been  filled  with 
data  earlier.  The  roles  of  the  two  buffers  are  reversed  when  the  real  time 
buffer  is  filled. 

3. 4 The  Fast-Fourier  Transform  (FFT)  Processor 

The  principles  of  operation  of  the  FFT  are  well  described  in  the 
literature.  Given  a time  series  input  N complex  (I  fc  Q)  samples  in  length, 
the  FFT  produces  an  N point  spectrum  analyzer  output. 

3.  5 Magnitude  Determination 

The  FFT  spectrum  output  is  complex,  comprised  of  I and  Q components. 
Only  the  magnitude  is  of  interest  to  us.  The  magnitude  unit  computes  a 
good  approximation  of  the  magnitude  J 1^  + for  each  spectrum  point 
(filter). 


Mi 
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3.6  Range-Cell  CFAR 

This  technique  is  well-described  in  the  literature.  See  References  (6) 
and  (7).  Each  FFT  filter  output  except  the  zero  frequency  one  has  its  own 
range-cell  CFAR  processor  which  produces  a threshold  which  applies  to  that 
filter  only. 

3. 7 Frequency-Bin  CFAR 

This  processor  sets  a threshold  for  each  range  cell  which  is  derived 
from  the  average  of  all  filter  outputs  except  those  around  zero  frequency.  As 
a result  any  signal  whose  energy  is  spread  more  or  less  uniformly  throughout 
the  spectrum  instead  of  being  concentrated  in  one  or  two  frequency  bins  is 
excluded.  This  gets  rid  of  certain  kinds  of  pulse  interference  and  jamming. 

3.  8 Clutter  Map 

Really  strong  point  clutter  can  break  through  the  MTl/FFT  process. 

The  clutter  map  suppresses  this  residue  by  storing  the  locations  of  these 
clutter  points  (which  appear  in  the  zero  frequency  filter)  and  setting  a 
threshold  in  the  zero  bin  and  proportionately  in  near-zero  bins. 

3. 9 Threshold 

We  have  described  three  sources  of  threshold  in  3.6,  3.7  and  3.8  above. 
The  signal  must  exceed  each  of  these  thresholds  if  it  is  to  be  output  to  the 
post  proces8or--cs8entially  a logical  "and". 
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SECTION  4.  ANALYTICAL  FTS  STUDY 


4.  1 Introduction 

This  section  reports  on  a clutter  map/CFAR  study  aimed  at  evaluating  the 
performance  of  several  candidate  false  target  suppression  (FTS)  signal  processors 
for  the  Experimental  Array  Radar  (EAR).  In  4.  2,  a theoretical  model  is  developed 
for  comparing  the  performance  of  several  digital  doppler  processors  in  realistic 
clutter  environments.  This  model  is  implemented  in  a computer  simulation  des- 
cribed in  Section  4.  3,  and  performance  curves  (probability  of  detection)  are  given 
in  Section  4.4  for  each  of  the  candidate  processors  and  different  signal-to-inter- 
ference  (S/l)  conditions.  From  these  curves,  the  advantages  of  including  the  MTI 
and  appropriate  pre-FFT  weighting  are  clear  for  heavy  clutter  environments. 

With  MTI  there  will,  however,  be  a tendency  to  lose  low-velocity  targets  or  those 
at  or  near  blind  velocities.  There  exist  techniques  (range-gated  MTI,  staggered 
prf)  which  can  minimize  these  effects. 

The  basic  components  of  the  signal  processor  are  given  in  Figure  4-^L.  - 
The  system  begins  by  sampling  the  inputs  from  the  analog  in  phase  (I)  and  quad- 
rature (Q)  channels  of  the  receiver  with  an  analog-to-digital  converter  (ADC)  of 
appropriate  characteristics.  The  resulting  output  is  fed  through  an  MTI  canceller 
ahead  of  the  FFT  in  order  to  reduce  dynamic  range  and  to  improve  clutter  rejection 
and  is  then  weighted  appropriately  to  reduce  doppler  sidelobe  effects.  Since  doppler 
processing  is  associated  with  operations  on  blocks  or  "batches"  of  samples,  the 
weighting  network  output  is  input  into  a buffer  storage  and  an  FFT  then  translates 
these  samples  from  the  time  domain  to  the  frequency  domain.  Finally,  the  FFT 
outputs  are  fed  into  an  FTS  processor  whose  function  is  to  reduce  the  number  of 
clutter  false  alarms.  This  is  achieved  by  applying  a number  of  threshold  tests 
to  the  range  gate/filter  amplitudes.  The  clutter  map  provides  a scan-to-scan 
history  of  the  clutter  in  the  zero  velocity  doppler  filter  and  defines  an  adaptive 
threshold  for  the  low  velocity  doppler  filters.  Also,  a cell-averaging  CFAR  processor 
further  eliminates  the  effects  of  clutter  in  the  higher  velocity  doppler  channels. 
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4.  2 Theoretical  Model 
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In  this  section,  a mathematical  model  is  presented  for  computing  detection  i 

probabilities  of  the  output  of  several  clutter  map/CFAR  (CM/CFAR)  signal  pro- 
cessors. Specifically,  the  configurations  analyzed  in  this  and  later  sections  are 

1.  ADC  -*  MTI  -»  W -*  FFT  * CM/CFAR 

2.  ADC  -*  MTI  .*  FFT  ^ CM/CFAR 

3.  ADC  -*  W .♦  FFT  .*  C^/CFAR 

4.  ADC  -*  FFT  CM/CFAR 

5.  ADC  -♦  MTI  -»  CFAR 

where  W denotes  a weighting  network.  Although  several  processors  are  con-  |i 

ij 

sidered,  one  general  math  model  is  derived  which  covers  each  of  the  specific  ! 

designs  with  appropriate  modification  of  the  signal  covariance  matrices  at  the  ! 

FFT  input.  , ■ 

4.  2.  1 Signal  Models 

Target,  ground  clutter,  weather  clutter,  and  noise  signal  models  needed 
for  evaluating  detection  performance  in  realistic  clutter  (Interference)  environ- 
ments are  given  below: 

a.  Target  Signal. 

The  target  model  appropriate  for  the  step- scanned  EAR  radar  is  the 
Swerling  I model.  Here  the  target  signal  power  returned  per  pulse  is  assumed  to 
be  constant  for  the  time  on  target  during  a single  scan,  but  to  fluctuate  indepen- 
dently from  scan-to-scan.  Expressed  in  statistical  terms,  the  normalized  auto- 
correlation function  of  target  cross  section  is  approximately  one  for  the  time  in 
which  the  beam  is  on  target  during  a single  scan  and  is  approximately  zero  for  a 
time  as  long  as  the  interval  between  scans. 

b.  Ground  Clutter  Model. 

As  described  in  Barton  [3],  ground  clutter  consists  of  a fixed  component 
due  to  returns  from  rigid  scattering  elements  within  the  clutter  patch  (e.  g.  build- 
ings,  fences,  power  poles,  etc.  ) plus  a fluctuating  component  due  to  returns  from 
moving  reflectors  such  as  swaying  trees,  vegetation,  etc.  The  clutter  voltage 
samples  for  a given  range  bin  are  therefore  in  general  characterized  temporarily 
by  a Rician  amplitude  distribution.  In  this  work,  however,  the  fixed  (DC)  component 
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of  clutter  is  assumed  absent.  This  assumption  is  reasonable  for  processors  con- 
taining an  MTI,  and  should  give  reliable  detection  probabilities  for  the  case  of 
point  clutter  and  no  MTI  in  all  but  the  first  few  doppler  bins.  Thus,  in  the  com- 
parisons between  systems  with  and  without  MTI,  the  detection  probabilities  for 
the  NO- MTI  case  may  be  somewhat  optimistic  for  the  first  few  FFT  filters  where 
the  given  detection  performance  should  be  considered  only  as  an  upper  bound  to 
the  actual  system  performance.  With  this  in  mind,  the  ground  clutter  voltage 
samples  for  a given  range  bin  are  assumed  to  be  characterized  temporally  by  a 
Rayleigh  amplitude  distribution  with  density  function 


-zZ/2o 


2 


(4-1) 


2 

Here  Pc  = 2 a is  the  average  clutter  power  from  a given  range  bin  and  is  propor- 
tional to  the  clutter  backscatter  coefficient,  a°. 


As  described  in  Boothe  [8],  since  the  clutter  backscatter  coefficient  varies 
from  clutter  patch  to  clutter  patch,  <j°  is  itself  a random  variable.  We  assume 
here  that  a°  is  a Weibull  variable  with  the  probability  density  function 


p(o°)  = j(o°)b_1  exp 
Cl 


(4-2) 


where 

b = — (A  - Weibull  slope  parameter  = 2) 


and  5°  is  the  backscatter  coefficient  averaged  over  the  CFAR  window.  In  reference 
[S'],  Boothe  shows  that  this  Weibull  assumption  is  reasonable  over  land  and  gives 
typical  values  of  the  slope  parameter  A.  For  semi-wooded,  hilly  terrain,  A = 2 
is  representative  of  the  average  clutter  environment  and  is  therefore  used  in  this 
analysis. 

Nathanson  [4]  shows  typical  land  clutter  profiles  versus  range  obtained 
oy  averaging  clutter  returns  from  many  pulses.  These  clutter  pro- 
files were  analyzed  by  Brooks  [9]  and  were  also  found  to  fit  reasonably  well  with 
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a Weibull  distribution  having  an  exponential  range  bin-to-range  bin  correlation. 
Equivalently,  clutter  power  versus  range  is  not  constant  but  exhibits  the  same 
spatial  variation  as  that  of  the  clutter  backacatter  coefficient.  Reference  [10] 
suggests  a strong  range  bin-to-range  bin  correlation  for  land  clutter  (i.  e. , p > . 99). 
In  this  analysis,  we  have  selected  the  correlation  p = l because  it  is  mathematically 
easy  to  evaluate  CFAR  detection  performance  in  this  case  and  the  performance  for 
p > . 99  should  be  very  nearly  identical  to  the  performance  with  p=  1.  This  assump- 
tion is  equivalent  to  assuming  that  the  average  clutter  plus  noise  power  is  equal  in 
each  range  bin  of  the  CFAR  window. 

Finally,  we  assume  the  pulse-to-pulse  correlation  between  ground  clutter 
returns  is  given  by  the  correlation  function 

- (2ttct  t)2/2  (I) 

R(t)  * e (4-3) 


where  oc  is  the  rms  frequency  spread  of  the  clutter  and  is  related  to  the  rms 
velocity  spread  of  the  scattering  elements  oy  by  the  familiar  doppler  equation 


At  C-Band  (X  = . 0545M),  the  spectral  width  of  ground  clutter  can  range  from 
about  . 6 Hz  on  a calm  day  to  about  18  Hz  for  winds  greater  than  40  knots. 

c.  Weather  Clutter  Model 

The  weather  clutter  is  assumed  to  have  a Rayleigh  temporal  distribution 
with  the  pulse-to-pulse  correlation 

R(t)  = exp  £ - (2tt  cwt)2/2  j coe(2iTfwT)  (4-5) 

where  a . analogous  to  a , is  the  spectral  width  of  the  rainstorm  which  at  C-Band 
ranges  from  about  67  Hz  (~ 6 ft/sec)  for  light  rain  to  about  145  Hz  (~  13  ft/ sec)  for 
heavy  rain  [3],  Also,  f denotes  the  center  frequency  of  the  storm  and  can  range 


(1) 


This  correlation  function  is  consistent  with  the  assumption  of  a Gaussian 
clutter  spectrum. 
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from  0 Hz  to  about  945  Hz  (~  50  knots).  It  should  be  pointed  out  that  the  correlation 
function  (4-5)  in  equivalent  to  the  assumption  of  a Gaussian  spectrum  for  the  rase 
of  zero  center  frequency,  f , and  although  it  deviates  from  this  assumption  Tor 
non-zero  center  frequency,  the  associated  spectrum  will  approximate  a Gaussian 
spectrum  for  frequencies  between  zero  and  l/2  PRF  (2.  5 kHz)  when  the  center  fre- 
quency is  large  compared  to  the  rms  frequency  spread  of  the  storm.  These  two 
extreme  cases  (center  frequencies  zero  and  945  Hz)  are  the  two  cases  analyzed  in 
Section  4.  3. 

Weather  clutter  will  generally  be  dispersed  in  range  over  several  cells.  In 
this  analysis,  the  mean  clutter  powers  over  the  CFAR  windows  to  the  left  and  to 
the  right  of  the  target  cell  are  assumed  constant  but  not  necessarily  equal. 

d.  Noise  Model 

The  noise  model  used  is  characterized  by  a Rayleigh  amplitude  distri- 
bution with  zero  pulse-to-pulse  correlation  (white-noise). 

4.  2.  2 Clutter  Map/CFAR  Concept 

The  basic  clutter  map/ CFAR  concept  used  here  for  a FTS  doppler  processor 
is  illustrated  in  Figure  4-2.  The  range  bin  of  interest  containing  the  target  is 
depicted  as  Rq  and  the  N range  cells  to  the  left  and  right  of  Rq  are  shown  as 
R-l’  • * ' R-N  and  Ri»  • • • RN’  respectively.  An  L-point  FFT  (filters  denoted  0 
th rough ■L'"l)  is  taken  of  the  signals  from  each  of  the  range  bins.  In  the  figure,  from 
the  FFT  on,  the  flow  is  shown  for  target  detection  in  the  kth  filter  (and  Rq  range 
bin)  only  but  a similar  flow  is  assumed  for  each  of  the  remaining  filters  (and  range 
bins)  except  for  the  zeroth  filter  where  only  a clutter  map  threshold  is  established. 
The  basic  idea  is  to  establish  appropriate  thresholds,  compare  the  kth  filter  output 
of  range  bin  Rq  with  the  largest  threshold  determined,  and  to  declare  detection  if 
the  signal  exceeds  this  maximum  threshold.  There  are  three  thresholds  estab- 
lished for  each  non-zero  filter  in  the  processors  studied  here,  two  CFAR  (constant) 
false  alarm  rate)  thresholds  and  a clutter  map  related  threshold.  (A  fourth 
threshold  derived  from  frequency  bin  averaging  - see  3.  7 - is  not  included  in  this 
analysis,  since  it  is  only  effective  under  conditions  of  pulse  interference  which 
are  not  included  in  the  input  signal  model.  ) 
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The  CFAR  thresholds  T^p^ , are  simply  weighted  averages*1*  of 

the  kth  filter  outputs  of  N range  bins  to  the  left  and  right  of  Rq,  respectively. 
These  two  thro  holds  are  intended  to  regulate  false  alarms  in  a noise  plus  weather 
environment.  The  leading  and  lagging  cells  are  averaged  separately  to  establish 
thresholds  that  will  regulate  false  alarms  in  a non. stationary  background,  as 
would  exist,  for  example,  when  the  window  confronts  the  edge  of  a weather  mass. 


A clutter  map  threshold  T^.^  is  established  for  the  zeroth  FFT  filter  (each 

range  bin)  to  regulate  ground  clutter  false  alarms,  and  represents  a weighted 

value  of  the  clutter  level  averaged  over  several  past  radar  scans.  A clutter  map- 
(K) 

related  threshold  T^.^  is  established  for  the  kth  filter  and  is  a fraction  of  the 
T°m  threshold.  Here  is  variable  and  is  programmed  into  the  processor 
according  to  the  expected  doppler  roll-off  of  the  ground  clutter. 


4.  2.  3 Detection  Probability- Model 

a.  Interference  Distribution  At  FFT  Output 

We  begin  by  establishing  the  noise  plus  ground  clutter  plus  weather 
clutter  signal  distribution  for  a given  range  bin  Rq  at  the  kth  FFT  filter  output 
Point  (B)  of  Figure  4-3.  To  do  this,  we  make  an  assumption  on  the  initial  distri- 
butions of  these  individual  interference  signals  at  the  output  of  the  A/D  converters  - 
Point  . Consistent  with  the  interference  models  of  Section  4.  2.  1,  each  of 

these  initial  signals  has  a Rayleigh  amplitude  distribution,  that  is,  the  associated 
in  phase  (I)  and  quadrature  phase  (Q)  signals  are  independent  Gaussian  variables 
with  zero  mean  and  appropriate  pulse-to-pulse  correlation,  the  degree  of  which 
depends  on  the  type  of  interference.  Accordingly,  we  can  develop  the  kth  FFT 
filter  interference  statistics  by  starting  with  a Gaussian  distribution  (representing 
any  one  of  the  interference  signals)  in  say,  the  I-channel,  following  this  distribution 
through  the  MTI,  weighting  and  FFT  networks,  and  finally  combining  the  resulting 
Gaussian  distribution  with  an  equivalent  Q-channel  distribution. 

Let  the  initial  I-channel  distribution  be  described  by  a zero-mean  Gaussian 
distribution  with  joint  density  - 


***  The  weighting  applied  is  consistent  with  a 10" ^ 
(PFA). 


probability  of  false  alarm 
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v1  = (v,,  v2 vn) 


is  the  transpose  of  an  nth  dimensional  vector  whose  elements  v^,  i = 1, . . . n denote 

I. channel  voltages  associated  with  the  first  n pulses,  respectively,  and  |MJ  is  the 

determinant  of  the  n x n covariance  matrix  describing  the  pulse-to-pulse 

correlation.  For  a sliding  window  3-pulse  MTI  and  the  inputs  v^,  . . . v^,  the  n-2 

outDUts  f - v - 2v.  + v.  i = 1,  . . . n-2  again  have  a zero  mean  normal  distri- 

*^  ii  i+li+2  _ 

bution  with  joint  density  function  given  by  (4-6)  with  v replaced  by  % and  Mn 

replaced  by  = AMfl  2A’\  where  Mn_2  is  the  matrix  Mr  with  the  1st  two 

rows  and  columns  deleted  and  A is  the  linear  transformation  matrix  defined  by 

l = A"v  (4-7) 


L n-2  x n-2 

If  a weighting  network  is  included  and  C.,  . . . C 7 denote  the  pre-FFT 

(2)1  n~C 

weights  (Taylor),  the  covariance  matrix  M _ for  the  zero  mean  I-channel 

n **  w 

Gaussian  distribution  at  the  FFT  input  is  given  by 

m(22  = cm(I)2c_1 

n-2  n-2 

where 


(4-9) 


For  convenience,  let  L = n-2  be  the  number  of  FFT  points,  W.t  . . . W denote  the 

(2\  1 1-1 
I-channel  FFT  inputs,  = [mijlLXL’  anc*  t*enote  t^ie  k**1  FFT  filter  output  in 

the  I-channel,  then 


Wy  * • 


(4-1 I) 


ly  - ^ W exp  [-i(2tr/L)(j-l)k]  , k=0, 

j » 

K(Ik)  -iJiy 

where  the  real  and  imaginary  parts  Rfy  and  Jiy  are  given  by 


my  - Y W.  cos  ^ (j-l)kj 


(4-12) 


J(y  - ^ W.  sin  (j-l)k 

. i L -J 


(4-13) 


Using  the  notation  N(^,  a ) to  denote  a normal  distribution  with  mean  (j  and  variance 
2 

n , it  is  clear  that 


my  * N(0,o‘k) 

•><y  e N(0,  ajk) 


(4-14) 


where 


2 

’Rk 

X 

m.  .C,  .C,  . 

lj  ki  kj 

i. 

1 

I, 

2 

’jR  ' 

l 

m. , d,  . d.  . 
ij  ki  kj 

j = 

1 

kj  = 

cos 

[xr 

Vi  " 

sin 

V,i-1,k] 

(4-15) 


In  a similar  fashion  for  the  Q-channel  distributions 


K(Qk)  e N(0,  o^k) 
•»(Qk)  e N(0.  0jk) 


(4-16) 


V 


Vv* 


where  again  and  o^k  are  given  by  (4-lb).  Now  combining  the  l and  Q channel 

• list ributions,  the  total  complex  FFT  output  at  the  kth  filter  becomes 

l\  L » iQ,. 


u(ik)  + j<Qkn  + i[R(Qk)  - jxyi 


(4-17) 


K(F  ) + i^(Fk) 

2 2 

whore  both  H ( F'^ ) and  J(F'k>  t N(0,  a^k  + an<*  are  independent.  Finally, 

assuming  a square  law  magnitude  unit  after  the  F FT the  combined  signal 

/.  R(Fk)2  + .»(Fk)2  (4-18) 


for  range  bin  j and  out  of  the  kth  FFT  filter  has  the  exponential  distribution 


p(Zjk)  = 2 exp 


(-  VV) 


where 


2 2 
aD1_  + ctuu 


(4-19) 


(4-20) 


Inns  much  as  the  distribution  (4-6)  is  representative  of  each  of  the  three  independent 
interference  sources,  the  distribution  (4-19)  is  also  representative  of  the  combina- 
tion of  ground  clutter,  weather  clutter,  and  noise  if 


2 2 2 
aNjk  + °gjk  + CTNjk 


(4-21) 


2 2 2 

Hero  the  noise,  ground  clutter,  and  weather  clutter  variances  o-,.,  ,a 

Njk  g]k  Njk 

respectively,  are  each  determined  from  (4-15)  and  (4-20)  after  the  computation  of 
the  appropriate  covariance  matrices  using  the  assumed  pulse-to-pulse  and  spatial 
correlations  given  in  Section  4.  2.  1 

4.  2.  4 False  Alarm  Probabilities 

For  a given  threshold  voltage  V,pk,  the  false  alarm  probability  at  the  kth 
filter  output  and  range  Rq  is  1 


'Although  a linear  magnitude  unit  is  actually  proposed,  a square  law  unit  is  used 
hero  to  simplify  the  mathematics.  Experience  shows  that  the  results  for  either 
will  bo  very  nearly  the  same. 


PFAk  = ~z  $"««P(-S/2pok2>d{ 


where 


exP  Ym,  * 


VTk/Zaok 


(4-22) 


(4-23) 


is  a normalized  threshold  equal  in  general  to  the  maximum  of  several  established 
values,  except  when  k = 0 and  it  is  simply  the  clutter  map  threshold.  In  this  case. 


(4-24) 


where  the  notation  ^ Z oq  is  used  to  denote  a sum  of  signals  from  range  bin  0 and 
filter  0 over  M-past  radar  scans,  and  Kj  is  a clutter  map  constant  associated 
with  a certain  false  alarm  probability.  Using  (4-19),  it  may  be  shown  that  the 
probability  density  for  Y.  is 


(Y,  1 / M \ M Y M-l  j MYbo  ] 

P,Y»°I  ^7,  (*;  ] Y»° 

Whence,  the  average  clutter  map  false  alarm  probability  is  given  by 


(4-25) 


PFA  = 
o 


1 + M 


(4-26) 


and  the  value  of  Kj  corresponding  to  an  average  false  alarm  probability  of  10  is 


K,  - M(106/M  - 1) 


(4-27) 


For  the  non-zero  filters,  a clutter  map  related  threshold  is  compared  with  the 
CFAR  thresholds  and  the  largest  is  chosen  for  the  detection  criterion.  Denoting 
the  normalized  clutter  map  related  threshold  by  A^,  we  set 
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Ak  " VTl/2ook  = Qk  VT</Zc 


(4-28) 


ft,  O 
k oo 


Ybo  = *k  YbO 


where 


ak  (°gok  /°goo  ) 


(4-29) 


is  the  expected  gound  clutter  roll-off  with  frequency.  The  probability  density  for 
A^  is  then  given  as: 


P<Ak'  = 


M-l  / MAj 

cxp  - kTsT 


and  the  average  false  alarm  probability  associated  with  A^  is 


<v  = 


(4-30) 


(4-31) 


With  regard  to  the  CFAR  thresholds,  we  let  and  denote  normalized 
thresholds  associated  with  the  range  bins  to  the  left  and  right  of  Rq,  respectively. 
Then, 


B*  IV 

ok  }=  -1 


(4-32) 


l h 


(4-33) 


ok  j = 1 


where  the  sums  are  over  the  N range  bins  to  the  left  and  right  of  Rq,  respectively, 
the  normalization  is  with  respect  to  the  average  interference  power  from  R^,  and 
k^  is  a CFAR  constant  associated  with  a given  false  alarm  probability.  It  is  now 
assumed  that  the  average  interference  powers  are  constant  over  the  cells  to  either 
side  of  Rq,  but  not  necessarily  equal  from  side  to  side.  Accordingly,  the  average 
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interference  power  on  each  cell  to  the  left  of  R it  denoted  2a  , , , and  to  the 

2 ® K 

right  of  RqI  the  notation  2a j ^ ia  uaed.  Paralleling  (4-30),  the  probability 
densities  for  and  are  given  by 


p<  ?k>  - 

1 fN)N( 

'L.' 

N-! 

( N 

t.k(N-l)J  ^kj  \ 

Jk- 

) 

exp  V ^"k 

where 

— 

f °-l,k2/oo,k2 

if 

«k  = 

Bk 

X 

II 

L °l.k2/ao,k2 

if 

Pk  = 

ck 

(4-34) 


(4-35) 


and  the  associated  false  alarm  probabilities  are 

1 


PFA. 


1 + Tkk2 
N 


TT 


(4-36) 


2 2 

Here  is  determined  assuming  o j ^ = °1  k ant*  *8  t*lere^ore  given  by 

k2  = N(.o-6/N-  A,,.2* 


(4-37) 


Finally,  on  the  average,  the  kth  filter  threshold  Y^k  is  given  as 


Ybk  = max(Ak’  Bk.  ck> 


and 


(4-38) 


•*« 

*'  -1 


PFA,  = 
k 


(i  ♦ "* 


'A’ 

M . 


M 


1 


TT 


if  Y 


bk 


otherwise 


(4-39) 


c.  Target  Detection  Probability  - Swerling  I Target 

2 

With  the  total  average  interference  power  2 oQ  k known  for  each  filter 
k,  the  average  signal-to-interference  powers  xQ,  * j,  . . . l vs  frequency 
out  of  the  magnitude  network  are  determined  by  first  passing  a sinusoid  of  unit  peak 
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amplitude  (representing  a normalized  fixed  amplitude  target  signal  over  the  L + 2 
pulses)  through  a 3-pulse  MT1  and  weighting  network,  if  used,  and  then  through  an 
L-point  FFT  and  finally  through  the  magnitude  unit.  Denoting  these  samples  by 
80* ' " ®L  l,the  avera8e  signal-to-interference  samples  vs.  frequency  for  a 
Swerling  1 target  are  given  as: 


">(£) 


(4-40) 


where  S/l  is  the  average  signal-to-interference  before  the  MTI; 
n2 

o,  k 


o,  k 

.2 


and  20  f is  the  total  interference  power  before  the  MTI. 


'1 


The  detection  probabilities  for  a single  pulse  out  of  the  magnitude  unit  are  then 
given  by  the  well  known  Swerling  formula 


l)k 


exp(-  Ybk/(1  + 


(4-41) 


and  the  average  detection  probabilities  are 


1 


. M 


Bkkl 

1 + 

M(  1 + x) 


if  Ybk  = Ak 


Dk 


(4-42) 


, N 


1 + 


Tkk2 
N(l  + x) 


otherwise 


As  remarked  earlier,  the  average  ground  clutter  powers  2a  . are,  for  each  k, 

gK 

constant  over  the  CFAR  window, but  are  nevertheless  random  and  proportional  to 
the  Weibull  distributed  backscatter  coefficient  a0.  Averaging  (4-42)  over  the 
Weibull  distribution  gives  the  overall  average  detection  probability 


<P 


Dk>  * 


(4-43) 


4.  3 Clutter  Map/C  FAR  Computer  Simulation 
4.  3.  1 Introduction 


A simulation  has  been  developed  for  computing  detection  probabilities  accor- 
ding to  the  model  given  in  Section  4.  2.  The  simulation  is  divided  into  two  separate 
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programs,  CFAR  and  CFARP.  A top  level  flow  diagram  and  listings  of  each  of 
these  programs  is  given  in  Appendix  A. 
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The  program  CFAR  generates  the  initial  interference  covariance  matrices  j 

and  performs  the  linear  transformations  necessary  for  establishing  the  mean  : j 

interference  powers  for  each  FFT  filter  at  the  output  of  the  magnitude  unit.  These 

i 

interference  signals  are  stored  on  direct  access  file  as  a function  of  frequency  i 

(filter  #)  and  are  automatically  called  for  by  the  program  CFARP,  which  in  turn 
computes  target  detection  probabilities  as  a function  of  target  doppler  frequency 
for  various  S/l  conditions. 

4.  3.  2 Simulation  Inputs 

The  operator  inputs  for  the  programs  CFAR  and  CFARP  are  given  in 
Tables  4-1  and  4-2  below.  j 


Table  4-1.  CFAR  Inputs 


Variable  Name 

Definition 

A 

Weibull  slope  parameter. 

AT  AY 

Sidelobe  level  control  parameter  for  Taylor  weights. 

FW 

Center  frequency  for  rain  storm  (Hz). 

GM 

Logical  variable  - Turns  ground  clutter  model  on  or  off. 
This  variable  should  always  be  set  = . True,  on  input. 

IWGT 

Logical  variable  - IWGT  = T for  Taylor  weighting; 
IWGT  = F for  no  pre-FFT  weighting. 

KCLUT 

Logical  variable  - KCLUT  = T results  in  a call  exit  since  a 
D.  C.  ground  clutter  component  has  not  yet  been  programmed 
into  the  simulation. 

MT1 

Logical  variable  - MTI  = T for  an  MTI,  and  MTI  = F for  no 
MTL 

MTIMDL 

MTIMDL  = 1 for  3 pulse  sliding  window  MTI;  MTIMDL  = 2 
for  narrow  notch  MTL 

MUWO 

2 

Backscatter  cross  section  of  the  rain  (M  ) - target  range  cell. 

MUW  1 

Backscatter  cross  section  of  the  rain  (M^)  - to  the  left  of 
target  cell. 

MUW2 

Backscatter  cross  section  of  the  rain  to  the  right  of  target  cell 

J 
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Variable  Name 


NBAR 


NM 

NPTS 

S1CSC 

SIGSW 

STGCDB 

STKCLDB 


STNDB 

THMDL 

WM 

WMDL 


Variable  Name 


1ATTN 

1DEALCF 

IFRESP 


\t 

i i 

IPROB 

t 1 

; 4 

f * 

ITATTN 

k ; 

E \ 

| 

iwgt 

tr 

MORCFAR 

Table  4-1.  CFAR  Inputa  (Continued) 
Definition 


Number  of  term*  used  in  Taylor  weighting  function  - used  only 
in  the  case  of  pre-  FFT  weighting 
Logical  variable  - Set  NM  = T in  all  runs. 

Number  of  FFT  points  (maximum  of  64). 

Spectral  width  of  ground  clutter  (Hz). 

Spectral  width  of  weather  clutter  (Hz). 

Signal-to- ground  clutter  ratio  (dB). 

Signal-to-point  clutter  ratio  (dB)  - not  used  in  present  form 
of  the  simulation  since  a point  clutter  model  has  not  yet  been 
added  to  the  program. 

Sign«rl-to- noise  ratio  (dB). 

Threshold  model.  Set  THMDL  = 1 in  all  computer  runs. 
Logical  variable.  WM  = T for  rain;  WM  = F for  no  rain. 
Weather  model  used;WMDL  = 1 in  all  computer  runs. 


Table  4-2.  CFARP  Inputs 
Definition 


Logical  variable.  IATTN  = T for  a plot  of  ground  clutter  plus 
noise  attenuation  vs.  frequency. 

Logical  variable.  IDEALCF  = T if  a plot  of  detection  proba- 
bility is  desired  for  a fixed  threshold. 

Logical  variable.  IFRESP  = T for  a plot  of  MTI  response  vs. 
frequency.  The  plots  are  normalized  to  the  noise  out  of  the 
Oth  filter. 

Logical  variable.  IPROB  = T for  Plots  on  Probability  Paper. 
Logical  variable.  ITATTN  = T for  a plot  of  clutter  map 
threshold  attenuation  vs.  frequency. 

See  CFAR  Inputs. 

Logical  variable.  Set  MORCFAR  = F if  only  the  ITATTN  or 
IATTN  plots  are  desired  as  output. 
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Table  4-2.  CFARP  Inputs  (Continued) 


Variable  Name 


Definition 


See  CFAR, 


MTINDL 

NCLMAP 


See  CFAR. 


NR CELLS 


PRED 


PRFA 


STIPL 


RHOG 


THMDL 


XMIN 


XMAX 


YMIN 


YMAX 


Age  of  clutter  map  in  number  of  interference  samples 
averaged. 

Total  number  of  range  cells  in  cell  averaging  CFAR 
NRCELLS/2  cells  to  each  side  of  target  cell). 

Logical  variable.  Set  PRED  = T if  predicted  clutter  map/ 
CFAR  detection  curves  are  desired. 

Logical  variable.  Set  PRFA  = T if  false  alarm  probability  at 
each  FFT  filter  is  desired  for  a clutter  map-related 
thresholding  system  only  (NO  CFAR). 

Logical  variable.  Set  STIPL  = T if  a plot  of  S/l  vs.  target 
doppler  frequency  is  desired 

Set  RHOG  = 1 in  all  runs.  This  implies  a unit  ground  clutter 
correlation  from  range  cell  to  range  cell. 

See  CFAR. 

Minimum  frequency  (x-coordinate)  for  plotting  routine 
Maximum  frequency  (x-coordinate)  for  plotting  routine 
Minimum  y-coordinate  for  plotting  routine. 

Maximum  y-coordinate  for  plotting  routine. 


CFARP 


s 


Y 


t 

i 

f 

i 

t, 


$INP1,  RHOG  = 1,  NRCELLS  = 32,  NCLMAP  = 6,  NPTS  = 64, 

MTIMDL  = 1,  ITATTN  = F,  IATTN  = F,  IFRESP  = F,  MORCFAR  = T, 

MTI  = T,  IWGT  = T,  STIPL  = F,  THMDL  = 1,  PR  FA  = F,  IDEALCF  = T, 
PRED  = T$ 

$INP2,  IPROB  = T,  XMIN  = 0.,  XMAX  = 5000.$ 

Finally,  in  specifying  varying  degrees  of  Taylor  weighting,  the  parameters 
ATAY  and  NBAR  defined  in  Table  4-1  are  required.  These  parameters  are  given 
below  in  Table  4-3  for  different  sidelobe  levels. 


Table  4-3,  Tavlor  Weight  Parameters 


Sidelobe  Level 

A 

NBAR 

-20 

. 9527 

3 

-30 

1.3197 

4 

-40 

1.  6864 

6 

-50 

2.  055 

9 

-60 

2.  42 

12 

-65 

2.  605 

14 

-70 

2.  785 

17 

-75 

2.  97 

19 

-80 

3.  15 

21 

-85 

3.  335 

22 

-90 

3.  52 

25 

4,  3,  3 Simulation  Outputs 
CFAR  Outputs 

The  CFAR  program  produces  very  little  printed  output  since  its  main  purpose 
is  to  generate  and  store  statistics  that  are  needed  by  the  CFARP  program.  Those 
outputs  printed  by  CFAR  are: 
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All  Namelist  Inputs 

SIGCSQ  ARRAY  (normalized^  ground  clutter  cross  section  each  FFT 
filter) 

SIGNSQ  ARRAY  (normalized  noise  power  each  FFT  filter) 

SIGWSQ  ARRAY  (normalized  rain  cross  section  each  FFT  filter) 

MUT,  MUCL,  KCL,  MUN,  MUCLF(2),  MU1,  MUW1=  MUW2 

where  MUT  = Target  cross  section  (M2),  MUCL  = total  ground  clutter  cross 
section,  KCL  = cross  section  of  point  clutter  (always  0),  MUN  = mean  noise  power, 
MUCLF  = cross  section  of  fluctuating  ground  clutter  component  (0)  and  the  re- 
maining variables  are  defined  in  Table  4-1. 

CFARP  Outputs 

The  main  CFARP  output  is  detection  probability  vs.  target  doppler  frequency 
for  given  signal-to-interference  inputs  and  a given  clutter  map/CFAR  signal  pro- 
cessor configuration.  Also  available  for  output  are  detection  probabilities  for 
fixed  threshold  (known  interference)  systems.  Such  a system  is  one  in  which  each 
FFT  filter  threshold  varies  from  range  bin  to  range  bin  but  is  constant  in  any 
given  bin  for  all  time.  A system  of  this  type  thus  defines  a limit  to  the  perfor- 
mance of  any  corresponding  CFAR  system  and  can  be  used  aB  a reference  for  es- 
tablishing a loss  in  S/N  due  to  CFAR  (CFAR  loss). 

Some  of  the  secondary  outputs  of  the  program  (which  were  included  primarily 
as  an  aid  in  the  initial  debugging  process)  include: 

e Clutter  map  false  alarm  constant 

e CFAR  false  alarm  constant 

e A plot  of  the  MTI  response  vs.  frequency.  The  plots  are  normalized  to 
the  noise  output  of  the  zeroth  FFT  filter 

• Plot  of  clutter  map  attenuation  vs.  frequency 

• Plot  of  S/l  vs.  target  doppler  frequency 


^ Normalization  is  with  respect  to  the  appropriate  cross  sections  at  the  signal 
processor  input. 

(2) 

The  MUCLF  variable  is  currently  not  used  by  the  program.  Since  there  is  yet 
no  point  clutter  model,  MUCL  represents  the  fluctuating  component.  Until  a 
point  plus  fluctuating  clutter  model  is  implemented,  a zero  will  be  printed  for 
MUCLF. 
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• Plot  of  ground  clutter  + noise  attenuation  vs.  frequency 
s Plot  of  false  alarm  probability  at  each  FFT  filter  for  clutter  map  related 
thresholding  (see  Section  4.  2.  3) 

4.  4 Simulation  Results 

4.  4.  1 Introduction 

In  this  section,  the  results  of  selected  computer  simulation  runs  are  pre- 
sented. The  objective  here  is  to  first  compare  the  detection  performance  of 
several  fixed  threshold  signal  processors  in  a heavy  ground  clutter  environment. 
Since  these  fixed  threshold  systems  define  a limit  to  the  performance  of  any  corres- 
ponding CF AR  system,  the  comparisons  can  be  used  to  define  the  basic  components 
and  parameters  of  a clutter  map/CFAR  signal  processor  that  would  be  effective  in 
comparable  clutter  conditions.  The  predicted  performance  of  such  a clutter  map/ 
CFAR  system  is  then  presented  and  compared  with  the  performance  of  the  corres- 
ponding fixed  threshold  system  for  various  S/N  and  S/C  conditions.  The  com- 
parison is  then  extended  to  include  the  effect  of  weather  clutter  on  detection 
performance. 

4.  4.  2 Performance  Curves 

Figure  4-4  shows  a comparison  between  each  of  the  fixed  threshold  systems: 

a.  No  MTI  and  no  pre-FFT  weighting 

b.  No  MTI,  but  40  dB  Taylor  weights 

c.  3 pulse  sliding  window  MTI  and  no  pre-FFT  weights 

d.  3 pulse  sliding  window  MTI  and  40  dB  Taylor  weights 

e.  Narrow  notch  MTI*  and  no  pre-FFT  weighting 

f.  Narrow  notch  MTI  and  40  dB  Taylor  weights 

g.  3 pulse  sliding  window  MTI  and  no  FFT. 

A 64-point  FFT  was  assumed  along  with  a Swerling  1 target,  S/N  and  S/C  condi- 
tions of  10  dB  and  -35  dB,  respectively,  and  ground  clutter  with  a spectral  width 
of  18  Hz.  Also,  the  curves  in  this  figure  reflect  an  average  false  alarm  probability 
of  10  ^ . It  is  seen  that  for  the  given  interference  conditions  a low  detection 


"Narrow- notch"  MTI  is  an  MTI  in  which  the  average  of  the  L pulses  required 
for  an  L-point  FFT  is  subtracted  from  each  pulse  before  entering  the  FFT.  For 
non-scanning  antenna  and  completely  non-moving  clutter,  this  type  of  MTI  is 
optimum. 
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probability,  reaching  a maximum  of  about  46%,  is  maintained  over  all  target 
doppler  frequencies  in  the  case  of  no  MTI  nor  pre-FFT  weighting.  Simply  adding 
40  dB  Taylor  weights  before  the  FFT  increases  the  detection  probability  on  the 
average  by  about  45  percentage  points  in  the  higher  filters,  and  a detection 
probability  of  more  than  80%  is  maintained  in  60%  of  the  filters  (target 
dopplers  between  1000  and  4000  Hz).  On  the  other  hand,  the  detection  prob 
ability  is  still  less  than  70%  for  target  dopplers  less  than  600  Hz  (i.  e. , for 
targets  appearing  within  the  first  8 FFT  filters). 

The  highest  detection  probabilities  shown  are  for  the  case  of  a 3 pulse  MTI 
and  no  pre  FFT  weighting.  In  this  case,  detection  probabilities  larger  than  97.  5% 
are  achieved  72%  of  the  time,  and  less  than  70%  detection  occurs  only  about  8%  of 
the  time  for  target  dopplers  less  than  about  200  Hz  (1st  3 filters). 

Following  the  MTI  with  40  dB  Taylor  weighting  degrades  detection  perfor- 
mance slightly  in  the  higher  filters,  but  improves  performance  significantly  in  the 
lower  filters.  Better  than  80%  detection  occurs  beyond  the  2nd  FFT  filter  and 
better  than  90%  detection  beyond  the  3rd.  The  degradation  in  detection  perfor- 
mance in  the  higher  filters  is  less  than  1%  and  is  expected  since  a reduction  in  S/N 
accompanies  pre  FFT  weighting. 

Also  shown  in  this  figure  is  the  significant  improvement  in  detection 
of  the  MTI/dopplor  processors  over  the  simple  MTI  systems  in  strong  clutter 
environment. 

In  order  to  interpret  the  results  of  Figure  4-4  in  terms  of  S/N  in  dB,  we 
resort  to  the  theoretical  developments  in  Section  4.  2.  From  Equations  (4-22)  and 
(4-41),  the  probabilities  of  false  alarm  and  detection  for  fixed  threshold  systems 
(Swerling  I target)  are  given  respectively  by  - 


PFA  = e 


- yb/(l  ♦ x_) 


1 

1 + xT 


(PFA) 


where  yb  is  the  threshold  and  Xp  is  the  signal-to- noise  ratio  required  in  the 
fixed  threshold  case.  Here  we  restrict  the  conversation  to  the  higher  FFT  filters 


where  ground  clutter  effects  are  negligible1 
the  expression  - 


For  these  filters,  (4-44)  leads  to 


general  Xp  is  the  average  S/L 
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3-Pulie  MT1:  No  Pre-FFT  Weighting 


Target  Doppler  Frequency  (Hz) 


Figure  4-4.  Probability  of  Detection:  Ideal  CFAR 
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(4-45) 


Ln  (PFA)  . 

XF  “ Ln  (PD)  ' 

for  the  signal-to- noise  ratio  required  for  given  false  alarm  rates  and  detection 
probabilities. 

Paralleling  the  above  for  CFAR  systems.  Equations  (4-26)  and  (4-42)  lead  to 
the  expression 


_ l 

- , (PFA)'.N_U.  ,4.46) 

C 1 

— “ N 
(PD)  1 - I 

N # Range  cells  averaged 

for  the  average  S/N  required  in  the  higher  filters  for  given  average  false  alarm 
rates  and  detection  probabilities.  Again,  it  should  be  pointed  out  that  these  ex- 
pressions  are  valid  for  Sterling  I targets  only  and  single  pulse  hits  out  of  the 
square  law  magnitude  unit  following  the  FFT. 

Equations  (4-45)  and  (4-46)  lead  to  the  curves  plotted  in  Figure  4-5  for  a 
false  alarm  probability  of  10"^.  Applying  the  fixed  threshold  curve  in  Figure  4-5  to 
the  curves  in  Figure  4-4  suggests  that  an  increase  of  more  than  14.  5 dB  in  signal* 
to- noise  is  achieved  in  the  highest  FFT  filters  by  adding  the  MTI  and  40  dB  weights 
to  the  no-MTI,  no-weighting  system.  Also,  the  degradation  in  detection  perfor- 
mance of  less  than  1%  (for  higher  filters)  when  40  dB  weighting  is  added  to  the 
MTI  system  corresponds  to  a loss  of  about  1.  2 dB  in  S/N. 

Also  indicated  in  Figure  4-4  are  results  for  a narrow  notch  MTI,  that  is, 
an  MTI  in  which  the  average  of  the  L-pulBes  required  for  an  L-point  FFT  is  sub- 
tracted from  each  pulse  before  entering  the  FFT.  It  is  not  surprising  that  there 
is  essentially  no  improvement  over  the  no-MTI  case,  as  the  notch  is  just  too 
narrow  to  eliminate  all  but  a D.  C.  component  of  ground  clutter.  It  also  does  an 
exceptionally  good  job  of  eliminating  the  targets  in  the  zero  doppler  bin.  Figure  4-6 
clearly  shows  this  to  be  true.  In  this  figure,  the  relative  responses  of  the  3-pulse 
and  narrow  notch  MTI  are  compared.  Also  shown  is  the  expected  roll-off  in  fre- 
quency of  18  Hz  ground  clutter. 
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Figure  4-5.  S/N  Requirements  for  Detection 


Relative  Response  (dB) 


Figure  4-6.  Relative  Response  Of  3-Pulse  and  Narrow-Notch  Filters 


Figure  4-7  shows  the  effect  of  changing  slightly  the  degree  of  pre-FFT 
weighting  for  the  fixed  threshold,  3-pulse  MTI  system.  By  varying 
the  40  dB  Taylor  weights  by  ± 10  dB,  the  signal- to- noise  requirements  for  the 
same  level  of  detection  in  the  higher  filters  vary  by  about  * . 5 dB.  That  is,  the 
CFAR  loss  is  ± . 5 dB.  Also  plotted  in  this  figure  is  the  detection  curve  for  a 
16  pt.  FFT  based  on  a single  pulse  out  of  the  magnitude  unit.  In  order  to  present 
a fair  comparison  between  the  16  and  64  pt.  FFT,  non-coherent  integration  of 
4 signal  plus  noise  variates  from  the  16  pt.  FFT  should  be  performed.  When  this 
is  done,  the  resulting  detection  curve  for  the  Swerling  1 target  is  raised  about  5.  5 
percentage  points  above  the  single  pulse  result  (see  Reference  [11],  page  238)  for 
the  higher  target  dopplers,  but  is  still  somewhat  below  the  detection  curve  for  the 
64  pt.  transform.  If,  however,  frequency  diversity  were  used  to  produce  rapid 
pulse-to-pulse  fluctuations  out  of  the  magnitude  unit,  Swerling  II  curves  from  [11] 
indicate  that  only  2 pulses  are  necessary  to  raise  the  detection  curve  for  the  16- 
point  transform  above  that  for  the  64  pt.  transform.  This  would  then  be  an  alter- 
native to  using  the  64-point  FFT  in  heavy  clutter  environments. 

Performance  curves  for  a clutter  map/CFAR  system  consisting  of  a 3-pulse 
MTI,  40  dB  Taylor  weights,  and  a 64-point  FFT  are  shown  in  Figure  4-8  along 
with  the  corresponding  fixed  threshold  case.  Also  shown  are  performance  curves 
for  simple  MTI  (No  FFT)  systems.  Here  the  same  heavy  clutter  environment  as 
in  thr  previous  figures  is  assumed.  For  doppler  frequencies  >50.0  Hz, 

the  CFAR  loss  in  dB  can  be  read  off  the  curves  in  Figure  4-5  . As  the  number  of 
range  cells  averaged,  each  side  of  the  target  cell,  decreases  from  32  to  8,  the 
CFAR  loss  increases  but  still  remains  less  than  1 dB.  Similar  performance  curves 
for  various  S/N  and  S/C  ratios  are  given  in  Figure  4-9. 

Finally,  Figures  4-10  and  4-11  show  clutter  map/CFAR  performance  for  an 
MTI,  40  dB  weighting,  64  pt.  FFT  system  operating  in  heavy  ground  clutter  and 
rain.  The  results  are  given  for  various  rainfall  rates  and  the  extreme  cases  of 
storms  with  center  frequencies  of  0 and  945  Hz. 

4.  4.  3 Remarks  and  Conclusions 

a.  In  heavy  clutter  environments,  a combination  of  MTI  and  appropriate 
pre-FFT  weighting  will  improve  the  detection  performance  over  a 
system  which  excludes  the  MTI. 
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Figure  4-9.  Clutter  Map/CFAR  Performance  Curves 
(Noise  k Ground  Clutter  Environment) 
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Figure  4-11. 


Clutter  Map/CFAR  Performance  Curves 
(Noise,  Ground  Clutter,  Weather  Environment) 


b.  Low-velocity  target*  or  target*  near  a blind  speed  are  lost  or  severely 
degraded  by  the  MTL  Several  thing*  can  be  done  to  improve  this 
situation. 

e Implement  a range -gated  MTI,  such  that  clutter  is  sensed  in  each 
range  cell  and  the  MTI  activated  only  if  the  clutter  amplitude  is 
significant. 

* Implement  a staggered  PRF  MTI  to  minimize  blind  speed  effects. 

• We  have  considered  two  MTI  configurations  here:  the  three-pulse 
MTI  and  the  narrow  notch  MTI.  The  three-pulse  MTI  is  much  better 
than  the  other  in  clutter  rejection  but  worse  in  the  rejection  of  low 
velocity  targets.  Neither  i*  really  optimum  for  the  assumed  target 
and  clutter  environment.  Further  study  is  needed  to  determine  the 
design  of  the  optimum  MTI  for  this  processor. 

c.  The  curves  presented  do  not  include  the  effect  of  strong  point  clutter  on 
target  detection.  Thought  should  be  given  toward  refining  the  simulation 
to  include  this  D.  C.  clutter  component. 

d.  The  simulation  effectively  handles  Swerling  I targets  based  on  one  hit 
only.  Non-coherent  integration  and  the  capability  of  handling  other 
Swerling  targets  should  be  included  as  future  refinements  to  the 
simulation. 

e.  At  present,  the  rain  model  is  reasonable  for  zero  center  frequency 
cases  and  cases  where  the  center  frequency  of  the  rain  storm  is  large 
compared  to  the  spectral  width  of  the  rain.  Some  thought  should  be 
given  toward  refining  the  model  to  include  intermediate  cases. 
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SECTION  5.  HARDWARE  CONFIGURATIONS  AND  TRADEOFFS 


5. 1 Introduction 

Having,  in  the  preceding  section,  determined  the  requirements  on  the 
FTS  processor,  we  are  left  with  the  task  of  translating  it  into  a hardware 
design.  There  is  a variety  of  choices  of  processor  architecture  available. 
This  section  discusses  the  choices,  their  hardware  requirements  (in  terms 
of  components),  software  requirements,  and  degree  of  development  necessary 
to  fit  them  to  the  EAR  needs.  All  of  the  approaches  discussed  here  can  meet 
the  system  requirements;  the  choice  is  then  to  be  made  on  the  basis  of 
relative  cost  and  design  flexibility. 

5.  2 Input  Interface 

The  A/D  converters  of  the  signal  processor  convert  the  bipolar  video 
radar  returns  (of  the  Inphase  and  Quadrature  channels)  to  a suitable  binary 
format  for  digital  processing.  These  converters  (existing  equipment)  are 
Computer  Lab  Model  HS-905  9-bit  converters.  They  will  be  operated  in  the 
EAR  at  the  rated  (sample)  5 MHz  rate  with  a dynamic  range  of  + 1024  MV. 

The  output  format  is  an  offset  binary  code  where 

1 1111  1111  denotes  + 1024  MV 

1 0000  0000  denotes  0000  MV 

0 0000  0000  denotes  - 1024  MV 

The  A/D  converter  outputs  and  the  input  code  commands  are  transferred 
to  and  from  the  Digital  Moving  Target  Indicator  (DMTI)*  logic  through  a 
digital  buffer  unit  consisting  of  digital  differential  line  drivers  and  receivers. 
The  output  from  the  A/D's  can  bypass  the  DMTI  (upon  command)  and  input 
directly  to  the  FTS  processor. 


This  refers  to  the  existing  DMTI,  rather  than  to  any  new  sliding  window 
DMTI  which  might  be  implemented  as  part  of  the  FTS  processor. 
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5.  3 Arithmetic  Precision  and  Dynamic  Range  Considerations 


$ 


The  existing  9-bit.  (8  bits  plus  sign)  A/D  converters  define  the  dynamic 
range  capability  of  the  processor;  allowing  2 or  3 bits  to  define  the  rms  noise 
level  leaves  5 or  6 bits  for  target  amplitude  variation,  implying  that  the 
maximum  signal-to-noise  ratio  that  the  processor  can  handle  linearly  is 
30  to  35  dB.  The  ratio  between  largest  linear  signal  and  smallest  detectable 
signal  can  be  greater  than  this,  however,  since  the  coherent  integration  due 
to  the  FFT  will  allow  signals  weaker  than  the  input  noise  level  to  be  detected. 
For  example  a 64-point  FFT  produces  18  dB  of  S/N  gain;  for  a 13  dB  output 
S/N  (consistent  with  reasonable  probability  of  detection  and  false  alarm  rate), 
the  input  S/N  is  -5  dB. 

We  can  expect  strong  point  clutter  to  sometimes  exceed  the  linear 
range  of  the  processor.  It  is  important  that  this  overload  be  handled  in  the 
right  way.  The  limiting  should  be  defined  in  the  IF  circuitry  of  the  radar, 
before  the  I and  Q phase  detectors.  If  the  A/D  converters  themselves  are 
allowed  to  be  overdriven,  undesirable  non-linear  effects  will  take  place 
(i.  e. , clutter  harmonics  being  generated  in  higher  frequency  doppler  bands) 
which  can  seriously  degrade  system  performance.  With  limiting  handled 
correctly  at  IF,  the  clutter  map  circuitry  in  the  FTS  processor  can  eliminate 
the  clutter  point  from  the  output. 

Although  the  processor  input  is  9 bits,  the  number  of  bits  within  the 
processor  grows  by  integration.  This  fact  plus  practical  considerations 
relating  to  the  bit  specifications  of  available  integrated  circuits  dictate  that 
internal  arithmetic  be  carried  out  to  16  bits,  although  FFT  coefficients 
and  sidelobe  weighting  coefficients  need  not  be  carried  higher  than  12  bits. 

5.  4 Input  Weighting 

The  weighting  function  is  multiplied  by  the  complex  input  data  stream 
before  it  is  stored  (as  a block  of  N samples  per  range  gate)  in  the  input 
buffer.  The  input  data  is  generated  as  an  offset  binary  code.  For  efficient 
multiplication  it  must  be  converted  to  a 2's  complement  representation.  The 
weighting  coefficients  are  stored  as  12-bit  numbers;  every  200  ns  two 
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10  x 12  bit  multiplications  are  performed  with  the  resulting  output  rounded 
to  a 10-bit  2's  complement  number  before  storage  in  the  buffer. 

5.  5 Input  Buffer 


i 


The  input  buffer  must  be  sized  to  handle  the  maximum  number  of  range 
bins  times  the  number  of  points  in  the  FFT.  Flexibility  requirements  dictate 
that  the  number  of  points  be  variable.  Then  we  will  need  to  vary  the  number 
of  range  cells  that  can  be  processed  as  well  if  the  buffer  size  is  to  be  constant. 
Each  range  cell  generates  a 20-bit  word  (10  bits  each  for  1 and  Q channels). 
Then  for  example  if  we  consider  768  range  gates  (bins)  the  maximum,  and 
the  minimum  transform  size  as  16  points/  then 
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) = 12,  288  words  x 20  bits 


We  choose  the  closest  binary  value,  16,384  words  x 20  bits  or  327,680  (this 
excess  will  provide  for  storage  of  intermediate  values  when  performing  the 
Fast  Fourier  Transform).  The  binary  approach  to  varying  the  range  gates 
also  keeps  the  control  mechanism  of  the  input  buffer  simple  when  writing 
into  or  reading  out  of  it. 

Since  the  PRF  is  5 kHz,  the  time  between  samples  is  200  usee.  For 
the  16  point  transform  the  memory  map  appears  as,  (with  the  time  to  collect 


Although  the  buffer  size  is  constant,  the  time  to  fill  up  a buffer  is  related 
to  the  number  of  points  in  the  transform  as  follows: 

768  range  gates  x 16  points  - 3. 2 ms 

384  range  gates  x 32  points  = 6.4  ms 

192  range  gates  x 64  points  = 12.8  ms 

96  range  gates  x 128  points  = 25.6  ms 

48  range  gates  x 512  points  = 102.4  ms 

24  range  gates  x 1024  points  = 204.  8 ms 

There  are  two  possible  techniques  which  can  be  used  for  the  input  buffer. 

5.  5.  1 Double  Buffer  Concept 

After  weighting,  the  digitized  1 & Q data  is  stored  in  one  buffer  while 
the  second  batch  of  data  is  fed  into  the  FFT.  This  is  depicted  below  and 
requires  twice  the  memory  size  or  2 x 16,  384  words  x 20  bits  = 655,  360  bits. 
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Buffer 
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The  advantage  of  this  technique  is  the  simplicity  of  control;  the  disadvantage 
is  that  of  requiring  twice  the  memory  size.  Further  detailed  examination  of 
the  read/write  speed  required  for  this  memory  and  its  implementation  is 
necessary.  The  Bipolar  memory  at  this  speed  (5  MHz)  comes  in  1024  x 1 size  and 
we  can  package  40  IC's  on  one  card,  or  about  2K  words  x 20  bit  on  a single 
8x10  inch  PC  board,  for  a total  of  16  cards  for  the  full  32K  double  buffer. 

Another  approach  is  to  multiplex  the  data  into  a high-speed  buffer 
(thus  widening  the  word  size)  and  transfer  the  wider  word  into  a slower 
speed  MOS  memory.  The  output  from  the  weighting  network  is  that  of  a 
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10-bit  number  every  100  n*.  If  we  distribute  this  by  multiplexing  into  two 
sets  of  ten  10-bit  registers  at  the  100  ns  rate  and  then  transferring  one  set 
into  a slower  MOS  memory  (100  bit  word  width),  we  can  perhaps  have  a less 
costly  solution  than  using  all  high-speed  memory. 

5.5.2  Single  Buffer  With  C orner -Turning  Memory 

This  technique  is  a single  buffer  concept  for  storing  the  batch  samples 
which  are  fed  to  the  FFT.  The  basic  corner-turning  idea  depends  on  the  fact 
that,  as  data  is  read  out  of  the  memory,  space  is  made  available  so  that 
other  data  can  be  written  in.  At  the  end  of  a read/write  sequence  however 
the  pattern  of  data  (i.  e. , the  address  locations  corresponding  to  a given 
range  ccll/sample  point)  has  changed.  The  new  pattern  is  however  known  to 
the  sequence  control  program,  and  the  whole  process  is  repeated  with  a 
new  set  of  addresses,  leading  to  still  a new  pattern.  An  so  on.  Hopefully, 
in  not  too  many  repetitions  the  original  pattern  will  repeat  itself. 

An  investigation  was  made  of  the  corner-turning  pattern  sequences 
involved  for  the  projected  memory.  It  was  soon  realized  that  the  requirement 
for  flexibility  made  the  addressing  problem  involved  in  corner -turning 
hopelessly  complex.  Although  for  certain  special  cases  simple  pattern 
repetitions  might  be  possible,  it  was  certainly  not  true  in  general.  A 
judgment  was  made  that  the  complexity  of  the  address  generating  circuitry 
would  be  sufficiently  high  in  a flexible  corner -turning  memory  so  as  to 
outweigh  the  saving  in  memory  involved  as  compared  to  the  double  buffer. 

5.  6 Fast-Fourier  Transform  (FFT)  Implementation 

5.6.1  General  Structure 

The  well-known  FFT  algorithm  efficiently  carries  out  the  Discrete 
Fourier  Transform  operation  described  by  the  following  equation: 


F(k)  = s (n)  e‘jnk(  n"  ) 


(5-1) 


whore 


! ? I 


f'-*K 


N 

s (n)  = 
k 

F (k)  = 


number  of  points 

complex  sample  corresponding  to  nth  point 
filter  number 

complex  spectrum  amplitude  of  kth  filter 


The  FFT  process  by  which  the  above  equation  is  partitioned  to 

minimize  the  necessary  computations  is  abundantly  described  in  the  literature. 

N 

Suffice  it  to  say  that  the  FFT  requires  about  — log2  N complex  arithmetic 
operations.  The  basic  element  of  the  FFT  is  the  fundamental  "Butterfly" 
operation  which  performs  the  necessary  complex  multiplications  and 
summations,  as  shown  below. 
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Figure  5-1.  Butterfly  Representation 


Note  from  the  above  equations  that  the  requirement  is  for  4 multiplication 
operations  and  6 add/subtract  operations.  The  four  multiplication  operations 
are: 


Ifi  x cos  0.  Qg  x cos  0,  Ig  x sin  0,  and  Qg  x sin0 


(5-3) 
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The  six  add/subtract  operations  are: 

Partial  Sum  #1  = Ig  cos  8 + Qg  sin  8 (5-< 

Partial  Sum  #2  = Qg  cos  8 - Ig  sin  0 

lA  + #1*  lA  ' #1*  °A  + #2'  QA  ' n 


Diagramatically,  the  above  mathematical  process  is  represented  as: 


Figure  5-2.  Physical  Representation  of  Butterfly 


Examining  the  above  diagram,  we  see  that  the  Butterfly  structure  can 
be  divided,  in  a hardware  organization  sense,  in  two  ways --vertical  slice 
or  horizontal  slice. 

Two  kinds  of  vertical  slice  organization  are  diagrammed  below.  There 
can  be  either  two  types  of  cards--Type  1 and  Type  2--or  all  Type  1 cards 
with  one  of  the  adders  not  used  in  half  of  the  cards. 


There  can  be  three  kinds  of  horizontal  slice  cards  as  shown  below. 
Or  we  can  get  along  with  only  Type  1 and  Type  3,  leaving  out  adders  when 
a Type  2 is  required. 
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[ffl 
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KK 


T ype  1 


f 


f 


Type  2 


Type  3 


5.6.2  Candidate  FFT  Approaches 


Three  general  approaches  to  the  FFT  hardware  design  will  be  considered 
and  compared.  Two  of  them  are  basically  existing  Raytheon  designs,  although 
they  must  be  modified  and  reprogrammed  to  fit  the  EAR  requirements.  These 
are  called  the  Universal  FFT  (UFFT)  and  the  General  Purpose  Signal 
Processor  (GPSP).  A third  design  is  built  around  the  use  of  common 


processing  elements  tied  together  by  a data  distribution  bus.  Although 
potentially  of  high  performance  and  flexibility,  this  approach  requires  more 
engineering  development  than  the  other  two. 

Still  a fourth  approach,  rather  specialized  and  somewhat  out  of  the 
main  stream  of  this  study,  involves  the  use  of  an  existing  Government-owned 
FFT.  The  characteristics  of  this  processor  and  the  problems  involved  in 
interfacing  it  with  the  EAR  are  discussed  in  Appendix  B. 

5.6.2. 1 The  Universal  FFT 

The  UFFT  design  {Raytheon  Equipment  Division)  is  representative 
of  the  vertical  slice  approach,  with  two  type  of  cards.  The  design  is  based 
on  a 2-bit  recirculation  (i.  e.,  2-bit  multiplier  and  2-bit  adder)  to  provide 
the  final  result.  The  clock  rate  is  8 MHz  (125  ns)  with  8 passes  through  the 
multiplier/adder  (8  x 125  ns)  to  provide  a complete  butterfly  in  1 Msec.  A 
recent  modification  to  this  design  has  enabled  a clock  rate  of  13.  3 MHz 
(73  ns),  speeding  the  butterfly  up  to  600  ns.  The  design  as  currently 
implemented  includes  the  scratch  pad  RAM  (for  immediate  storage  of  the 
two-bit  partial  results)  and  the  control  cards  for  a total  of  12  cards 
(8"  x 6"). 

5.6.  2.  2 The  General  Purpose  Signal  Processor 

Another  approach  to  the  FFT  implementation  is  the  GPSP  design 
(Raytheon  Missile  Systems  Division).  This  implementation  provides  a 
butterfly  in  200  ns,  using  a pipeline  construction.  This  architecture,  shown 
in  Figure  5-3,  is  a modification  of  the  horizontal  slice  approach.  This 
approach  requires  a minimum  of  50  cards. 

5.6.2.  3 The  Common  Element  Approach 

Another  approach  is  to  develop  a common  processing  element  (a  micro- 
processor) which  can  do  all  of  the  necessary  FFT  operations  for  a certain 
number  of  range  cells.  As  more  range  cells  are  needed,  more  common 
elements  are  added  to  the  bus -interconnected  system.  This  approach  has  a 
high  degree  of  modularity  and  flexibility.  It  is  diagrammed  in  Figure  5-4. 


i 


OUTPUT 


OUTPUT 


Figure  5-3.  Programmable  (Pipeline)  Signal  ProceBBor 
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The  rationale  for  this  scheme  derives  from  the  fact  that  in  the  FFT 
the  I/O  time  is  small  compared  to  computation  time.  The  common  element 
approach  also  has  this  behavior.  The  processing  time  per  range  gate  versus 
the  number  of  points  in  the  transform  is  given  by  the  equation, 

processing  time  = nun?-fr?r  S£-P°^nts  fo)  x log2n  x time  per  Butterfly  (tfip) 

The  I/O  time  as  a function  of  the  number  of  points  (n)  is  given  by  the 
equation, 

I/O  time  - n x 2 x time  for  I/O  cycle  Uj/q)  (5-5) 

A third  quantity  of  interest  is  the  number  of  common  elements 
necessary  to  perform  the  FFT  process  (as  related  to  the  Butterfly  time) 


as  a function  of  the  number  of  points,  and  is  found  from  the  following: 


number  of  range  gates/element  = -E^-8e  repetition  interval  (time)  x n 

processing  time  for  a single  element  ' ' 


number  of  elements 


number  of  range  gates  in  system 
number  oi:  range  gates/element 


(5-7) 


A plot  of  processing  time  versus  the  number  of  points  in  the  transform 
with  tfiF  as  a parameter  is  shown  in  Figure  5-5.  A plot  of  I/O  time  versus 
the  transform  size  is  shown  in  Figure  5-6.  The  two  plots  used  together  can 
show  the  ratio  of  I/O  time  to  processing  time  for  the  FFT.  For  a 4 Msec 
Butterfly  time  and  a 60  ns  cycle  time  of  the  memory: 
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1 


i 


No.  of  pt s I/O  Time 


Processing  Time 


I/O  Time 


Processing  Time 


16 

1.92  Us 

128  us 

1.  5% 

32 

3.  84  us 

320  us 

1.2% 

64 

7. 64  us 

768  us 

1.0% 

128 

15.3  Us 

1792  us 

0.  85% 

256 

30.  72  us 

4096  us 

0.75% 

512 

6 1. 44  us 

9216  us 

0.  67% 

1024 

122. 88  us 

20480  us 

0.6% 

Using  the  common  element  approach  to  perform  the  FFT  process  then  for 
768  range  gates,  a PRF  of  5 kHz,  and  a 16-point  transform  then: 


. . 768  range  gates  

number  of  elements  = 25  ^n^gates/Sl'eme*  = 30*  72 


31  elements 


Using  high-speed  microprocessor  techniques,  one  element  can  occupy 
one  card.  This  approach  if  it  can  be  properly  realized  in  hardware,  can 
represent  a significant  saving  in  total  hardware. 


5.7 


Signal  Processing  Hardware  Requirements 

This  section  summarizes  the  hardware  requirements  as  a function  of 
(a)  The  processing  approach  (GPSP,  Common  Element,  UFFT/C PE*),  (b) 

the  number  of  range  gates  processed,  and  (c)  the  number  of  FFT  points 
processed.  The  various  functions  which  need  to  be  performed  by  the  EAR 
Processor  are: 


i 

r 

1 

1 

(a) 

Input  Buffering 

(e) 

Range  Averaging  CFAR 

+ 

t 

(b) 

Weighting 

(f) 

Frequency  Averaging  CFAR 

* 

(c) 

FFT 

(g) 

Thresholding 

(d) 

Magnitude  Determination 

(h) 

Clutter  Mapping 

number  of  range  gates/element  = ^ ' = 25  ranges  gates/element  (5-8) 

1 b U |m|  S vC 


(5-9) 


The  CPE,  or  Common  Processing  Element,  approach  is  a microprocessor- 
based  technique  developed  at  the  Equipment  Division  for  handling  CFAR-type 
functions.  It  has  been  used  in  conjunction  with  the  UFFT. 
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The  GPSP  and  CE  (Common  Element)  approaches  can  do  all  of  these 
except  (a)  and  (h).  The  UFFT  does  (b),  (c)  and  (d).  The  CPE  does  (e),  (f) 
and  (g).  (a)  and  (h)  are  basically  memory  card  requirements  with  some 

processing  (of  the  CE  or  CPE  type)  also  needed  in  (h). 

The  way  in  which  the  various  functions  are  handled  by  each  of  these 
approaches  is  illustrated  in  the  functional  block  diagrams  of  Figures  5-7, 
5-8,  and  5-9.  The  functions  are  not  explicitly  shown  in  the  Common  Element 
Diagram  (Figure  5-9)  since  all  elements  are  basically  the  same,  and  are 
reprogrammed  as  needed  to  fulfill  the  necessary  functions. 


For  each  approach,  relationships  exist  between  the  number  of  cards 
needed  for  each  function  and  such  system  parameters  as  number  of  FFT 
points  and  number  of  range  cells.  These  relationships  are  Bet  down  below; 
they  are  derived  from  known  and  estimated  characteristics  of  these 
processors.  Memory  and  clutter  map  requirements  are  common  to  all 
approaches  and  are  discussed  separately  later. 


Symbols 

N 

R 

I 

T 

o 

B 

C 

T 

P 

INT 


number  of  FFT  points  processed 
number  of  range  cells 
number  of  operations/second 
scan  time 

number  of  FFT  butterflies 

number  of  cards  (40-60  IC's  each) 

-4 

processing  time  = 2x10  N seconds  for  a 5 kHz  PRF 
stands  for  "integral  part  of' 


5.7.1  The  Common  Element  Approach 

5.7. 1.1  Weighting 

Requires  operations  as  follows;  (per  range  gate) 


3 

1 

1 

I 

3 


j 

i 
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Figure  5-8.  Signal  Processor:  GPSP  Architecture 


BUS 


Distributed  Signal  Processor 
(Using  Common  Elements) 


Weighting  Process 

LOAD  Q term  in  Latch 

MPY  Weighting 

STO  Q term  in  place 

LOAD  I term  in  Latch 

MPY  Weighting 

STO  I term  in  place 

. Insert  Butterfly  process 

MPY  Weighting  (2nd  term  of  Butterfly  - I) 


MPY  Weighting  (2nd  term  of  Butterfly  - Q) 


Operations 


Index  Weighting  Address 


T c*al 


53  operations 


The  total  time  for  a block  of  N samples  at  a 5 kHz  PRF  is 
Tp  = 2 x 10_4N  seconds.  Then  for  R range  gates,  the  operations/sec.  1^ 
required  for  weighting  is: 


53R/T 


5.7.  1.2  FFT  Processing 


265,000  R/N 


(5-'  10) 


The  total  number  of  Butterflies  B per  range  cell  is: 

B = y-log2N=  1. 66  NlogjgN 

The  operations  per  Butterfly  are  tabulated  below: 
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FFT  Process  (Operations) 

Using  a 16  Sample  Batch  as  an  Example 


Operations 


LOAD 

Ig  in  Latch 

Latch 

1 

MPY 

Cos  8 

lg  cos  8 in  Accum 

12 

STO 

Temp  1 

Temp  1 = Ig  cos  8 

1 

MPY 

Sin  8 

Ig  cos  8 in  Accum 

12 

STO 

Temp  2 

Temp  2 = ' ^ sin  8 

1 

LOAD 

Qg  in  Latch 

Latch 

1 

MPY 

Cos  8 

Qg  cos  8 in  Accum 

12 

STO 

Temp  3 

Temp  3 = Qg  cos  8 

1 

MPY 

Sin  8 

Qg  sin  8 in  Accum 

12 

ADD 

T emp  1 

Accum  = Ig  cos  8 + Qg  sin  8 

1 

STO 

Temp  1 

Temp  1 = Ig  cos  8 + Qg  sin  8 

1 

LOAD 

Temp  3 

Latch 

1 

SUB 

Temp  2 

Accum  = Qg  cos  8 - Ig  sin  8 

1 

STO 

Temp  2 

Temp  2 = Qg  cos  8 - Ig  sin  8 

1 

LOAD 

*0 

Latch 

1 

ADD 

Temp  1 

Iq  = Iq  + (Ig  cos  8 + Qg  sin  8) 

1 

STO 

*0 

Store  in  original  location 

1 

SUB 

Temp  1 

Iq  now  in  Accum 

1 

SUB 

Temp  1 

Accum  = Iq  - ^Ig  cos  8 + Qg  sin  8) 

1 

STO 

X8 

Store  in  original  location  now  Ig 

1 

LOAD 

Q0 

Qq  in  Latch 

1 

ADD 

Temp  2 

Accum  = Q0  + (Q  cos  8 - I sin  8) 

1 

STO 

Q0 

Store  in  original  location  now  Qq 

1 

SUB 

Temp  2 

Qq  in  Latch 

1 

SUB 

Temp  2 

Accum  = Qq  - (Qg  cos  8 - Ig  sin  8) 

1 

STO 

Q8 

Store  in  original  location  now  Qg 

1 

Index 

*0 

1 

Index 

Q0 

1 

Index 

*8 

1 

Index 

°8 

1 

Total  74  operations 


As  above,  we  find  the  total  operations  as  follows: 

■f  ■ 6.14  x 10^R  logj0N 

5.7. 1.3  Magnitude  Processing  (per  Range  Gate) 


The  magnitude  calculation  is  based  on  a 5 -level  approximation  algorithm 
which  performs  the  magnitude  calculation  based  on  the  values  of  I fit  Q from  the 
point  of  the  output  of  the  FFT  process.  This  algorithm  produces  an  output  with 
a peak  error  no  worse  at  any  point  than  3.  1%  of  the  correct  J 1^  + value. 
Both  a table  of  operations  and  a flow  chart  of  the  magnitude  processing 
algorithm  are  shown  below. 


i 

i 


Magnitude 

Process 

Operations 

LOAD 

I 

I in  Accum 

1 

SUB 

Q 

I-Q  in  Accum 

1 

BMI 

SWAP 

I<  Q? 

1 

LOAD 

I 

1 

STO 

Temp  1 

store  larger 

1 

LOAD 

Q 

1 

LARG: 

STO 

Temp  2 

store  smaller 

1 

SHR 

shift  right 

1 

SHR 

divide  by  4 

1 

ADD 

Temp  1 

1 

STO 

Temp  3 

1 

LOAD 

Temp  1 

1 

SHR 

1 

STO 

Temp  1 

C/2 

1 

SHR 

1 

ADD 

Temp  1 

1 

STO 

Temp  4 

3/4  C 

1 

SUB 

Temp  2 

3/4  C-D 

1 

BMI 

EIT 

D > 3/4  C 

1 

LOAD 

Temp  1 

C/2 

1 

SUB 

Temp  2 

C/2-D 

1 

BMI 

THTO 

D >C/2 

1 

LOAD 

Temp  3 

Load  magnitude 

1 

FINAL: 

STO 

in  I 

store  in  I 

1 

EXIT 

1 

SWAP: 

LOAD 

Q 

Q in  Accum 

1 

STO 

Temp  1 

store  larger 

1 

LOAD 

I 

load  smaller 

1 

BR 

LARG 

return 

1 

EIT: 

LOAD 

Temp  2 

load  smaller  (D) 

1 

SHR 

1 

SHR 

1 

SHR 

divide  by  8 

1 

ADD 

Temp  3 

mag  + D/8 

1 

BR 

FINAL 

1 

THTO: 

LOAD 

Temp  2 

load  smaller,  D 

1 

SHR 

1 

SHR 

1 

SHR 

1 

SHR 

1 

SHR 

i 

divide  by  32 

1 

ADD 

Temp  3 

mag  + D/32 

1 

BR 

FINAL 

1 

EXIT: 

Index  I 

1 

Index  Q 

1 

Worst  Case  Total  34  operations 
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The  magnitude  operations  per  second  are  independent  of  block  size, 
and  is  given  by: 

L.  = 34  NR/(T  N)  = 170,  000  R (5-13) 

M p 

5.  7.1.4  Range  Averaging  CFAR 

The  Range  Average  processing  is  based  upon  summing  16  range  cells 
before  the  range  cell  of  interest  and  16  summing  range  cells  after  the  range 
cell  of  interest.  A range  cell  space  is  inserted  on  both  sides  of  the  range 
cell  of  interest  and  the  two  sums.  A tabulation  of  operations  follows: 


Range  Average  Process  Operations 


BR 

No  op 

; if  count  less  than  16 

1 

LOAD 

s i 
n- 1 

1 

SUB 

R 

n 

1 

ADD 

Rn+ 1 

1 

STO 

S 

n 

; late  sum 

1 

CMP 

S„-19 

; compare  both  sums 

1 

BR 

LARG 

1 

LOAD 

SUM 

; load  larger  sum 

I 

SHR 

1 

SHR 

1 

SHR 

1 

SHR 

; divide  by  16 

1 

STO  THRESHOLD 

; store  in  threshold 

1 

INDEX 

S 

n 

1 

INDEX 

Sn-19 

1 

INDEX 

Rn+ 16 

1 

INDEX 

COUNT 

; total  number  of  range  gates 

1 

BR 

COUNT  GREATER 

1 

Total 

20  operations 

Ope  ration 

rate  here  is 

given  by: 

Ira  = 1000,000  R 
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5.7.  1.5  Frequency  Averaging  CFAR 

This  is  basically  a broad-band  interference  reduction  technique.  It 
is  accomplished  by  summing  the  magnitude  from  all  the  filters  of  the  output 
of  the  FFT  Process  and  dividing  by  the  total  number  summed.  Tabulating 
ope  rations: 


Frequency  Average 
(tor  a 64  point  transform) 
LOAD  f 

ADD  fj 


Operations 


1 

1 


• 

• 

ADD 

f63 

Accumulate  64  terms 

1 

SHR 

SUM 

* 

1 

SHR 

SUM 

1 

SHR 

SUM 

1 

SHR 

SUM 

1 

SHR 

SUM 

divide  by  64 

1 

SHR 

SUM  ' 

1 

CMP 

THR 

compare  to  threshold 

1 

STO 

LARGER 

1 

Total  72  operations 


Operation  rate  1^,^  = 6875R 
5.7. 1.6  Thresholding 

This  operation  is  based  upon  the  following  computation:  four 

thresholds  are  compared  (clutter,  range  average,  frequency  average  and 
minimum  threshold),  the  largest  is  selected  and  multiplied  by  a CFAR 
constant,  then  compared  to  the  magnitude  of  the  range  gate  of  interest. 


Threshold  Process 

LOAD  Clutter 

CMP  Range 

BR  Larger 

CMP  Freq. 

BR  Larger 

CMP  Min. 

BR  Larger 

MPY  CFAR  Constant 

CMP  Target 

FORMATTING 


I.I  H 200,000  H/N 


Operations 

1 

1 

1 

1 

1 

1 

1 

12 

1 

20 

Total  40  operations 


I 
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The  total  of  all  of  the  above  operation  rates  is  given  by: 

LpoT  = ^C  276'875  + 4-6^°°-  + 6.  14  x 105  log  (! 

Given  a card  implemented  with  high-speed  micrologic  which  is  capable 
of  processing  16.67  MIP  (16.67  x 10^  = I),  the  number  of  cards  required 
for  all  operations  is: 


JCE 


INT 


( il ) 

\ 16.67  x 106  / 


+ 1 


(! 


Since  one  card  is  needed  for  control,  etc. 

5.7.2  The  GPSP  Approach 

Since  this  is  an  existing  design  whose  properties  are  known,  we  do 
not  need  to  approach  it  analytically  as  above.  We  simply  note  that  50  cards 
as  a unit  can  process  10  butterflies/sec.  , as  well  as  handling  the  other 
needed  operations.  But  the  required  number  of  butterflies/ sec.  is: 


B 


KEU 


B = 8.  3 x 103  R log  N 


(! 


! ! 

? 

i 

i 

t 


The  GPSP  can  only  be  used  in  modules  of  50  cards.  So  the  card  count  here 
is  given  by: 


C(  i PS  P 


™ 0(^1  * •] 

50  Q INT  ( 8.  3 x 10*4  R log  N)  + 1 [] 


(! 


5.7.  1 The  UFFT/CPE  Approach 


For  the  UFFT,  12  cards  process  1.67  x 10  butterflies/sec.,  so  as 


above: 


c; 


UFFT 


12 


£lNT^4.  97  x 10*3  R log  + ij 


+ 2 


(! 


(The  2 represents  array  multiplier  cards  needed  for  weighting.  ) 
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-14) 


-15) 


-16) 


-17) 


-18) 


The  CPE  total  instructions  per  second  are: 

I-ICPK  *RA  + *FA  4 ^ H (calculated  as  in  5.7.  1) 


R £ 106,  875  + 20Pi  00°.  J 


Since  the  CPE  processes  3.  3 x 10°  operations/sec.,  and  requires  about  1 
control  card  lor  each  six  processing  cards: 

‘ err  h ) * m * 2 

\3.3X106/  \ 2X107/ 

The  total  card  count  for  the  UFFT/CPE  approach  is  the  sum  of 
CUFKT  and  CCPE* 

6 .7.4  The  Input  Buffer 

The  memory  requirement  is  given  by: 

tn  R x N x 2 (for  I & Q)  x 2 (for  double  buffering) 

- 4RN  words 

Assume  that  the  most  that  we  can  put  on  one  card  is  2048  words 
( 16- bit  words,  1024  x 1 memory  IC's,  32  memory  IC's  per  card  plus  some 
control  circuitry).  Then  the  number  of  buffer  cards  needed  is: 

‘‘•IB  ‘ INT  (#£)  ♦ 1 

6.7.5  The  Clutter  Map 

Assume  50  angle  cells  (100°  scan  in  2°  steps)  and  R range  gates, 
with  16  bits  of  averaged  clutter  amplitude  information  stored  per  cell. 

(This  is  probably  too  much,  and  will  lead  to  a conservative  estimate.  ) Then 
with  2048  words  stored  per  card,  the  number  of  cards  needed  for  clutter 
map  memory  alone  is: 

CCM  = mT  ( life  ) * 1 
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The  clutter  map  calculation,  for  a decaying  exponential  clutter 
memory,  will  follow  this  algorithm, 

where  K = 1-2 


S 


n+1 


K S + 
n 


lRil 


Where  R.  is  the  clutter  amplitude  of  a newly  arriving  sample,  SR  is  the 
stored  clutter  sum,  and  K is  an  exponential  weighting  constant.  The 
operations  are  tabulated  below: 


Clutter  Map  Process 


LOAD 

SHR 

SHR 

SHR 

SHR 

SHR 

SHR 

SUB 

ADD 

ST  O 

INDEX 


n 


R. 

1 


n 


old  value 
divide  by  64 


old  value 

add  magnitude  of  range  cell 
new  value  of  Sn 
update  for  new  sum 

Total 


Operations 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


11  operations 


The  map  processing  load  is  a function  of  the  scan  time  T^  (nominally 
0.  S second). 


lCMPROC  ' 50  angle  cells  x R 

- 1100R  ops/sec. 


x 11  /0.  5 sec. 

range  cells  ops 


Assuming  that  the  CPE  approach  (3.  3 MIPS)  is  used  for  implementing  this 
algorithm,  we  have  approximately: 

' 1100R  \ 

3.  3 x 106 


'CMPROC 


INT 


+ 1 


S.7.6  Other- -The  Synchronizer 

Six  (6)  cards  are  assumed  for  the  synchronizer,  common  to  all 
approaches. 


-23) 


6-24) 


(6-25) 
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5.7.7  Hardware  Comparisons 


A comparison  of  hardware  cost  alone  (and  we  can  consider  this  as 
being  proportional  to  total  card  count)  is  of  interest  in  that  it  is  a measure  of 
the  relative  cost  of  signal  processors  beyond  the  first  model  (within  which 
are  concentrated  development  and  software  costs).  Figure  5-11  gives  an 
exam  pie  of  the  relative  "cost"  measured  in  cards,  of  a processor  with 
N - 64,  for  varying  range  cell  number,  and  for  the  three  processor 
approaches  considered  here.  The  Common  Element  approach  appears 
consistently  to  cost  less,  although  not  always  by  a significant  amount. 

These  curves  have  been  plotted  from  the  equations  presented  in  the 
preceding  paragraphs.  Another  useful  way  of  looking  at  these  results 
appears  in  Figure  5-12;  here  we  consider  the  situation  of  a constant  cost 
system  (C  - 100),  and  examine  what  performance  is  available  (measured 
by  range  gates  R and  points  N)  using  the  three  approaches.  Again  the 
CE  approach  has  a slight  edge  except  where  the  number  of  range  gates 
becomes  of  the  order  of  1000  or  more,  where  the  GPSP  appears  to  gain 
an  advantage. 

5.  8 Development  Requirements 
5.  8.  1 Software 

The  various  operation  sequences  tabulated  in  the  preceding  paragraphs 
do  more  than  place  a demand  on  the  hardware  in  terms  of  operations  per 
second.  They  are  programs,  and  must  be  written  and  debugged  for  all  of 
the  approaches  considered.  Summarizing  the  lines  of  code  that  need  to  be 
programmed: 


Input  Weighting 
FFT  Operation 
Magnitude 
Frequency  Average 
Clutter  Threshold 
Range  Average 
Thresholding 

Multiply  Subroutine 


50%  Margin 


Total 


9 

30 

45 

72 

11 

20 

40 

227 

14 

241 

120 

36 1 lines 

x 4 each  line  equivalent  to  four 
1444  16 -bit  lines 
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Assuming  that  each  fully  documented  and  debugged  line  represents  the 
expenditure  of  1 engineering  man-hour,  the  total  software  cost  for  all 
approaches  is  about  200  man-days. 

5.8.2  Other  Costs 

(a)  Hardware  Development:  Even  when  card  designs  exist,  they 

must  be  packaged,  interconnected,  furnished  with  power  supplies  and 
controls,  etc.  For  the  Common  Element  (CE)  approach,  the  card  designs 
are  new  and  must  also  be  accomplished. 

(b)  System  Design  and  Test:  This  includes  working  out  many 

system  design  details,  writing  appropriate  unit  and  interface  specifications, 
and  testing  the  finished  unit  (which  usually  also  implies  the  design  of 
special  testing  hardware).  This  also  includes  documentation  of  the 
finished  unit. 

The  hardware  cost  data  of  the  preceding  section  was  developed  in 
terms  of  equivalent  "cards".  Most  of  the  development  fixed  costs  are 
basically  engineering  manpower  costs.  So  as  to  be  able  to  use  the  original 
curves  without  modification,  conversion  factors  are  used  to  relate  "cards" 
to  man-days.  Using  the  neutral  "card"  unit  seems  permissible  because 
relative  costs  are  of  prime  interest  here.  Also,  using  dollars  can  be 
misleading  because  of  inflationary  factors.  The  basic  equivalence  is  1 
tested  card  = 8 engineering  man-days. 


Summarizing  the  fixed  development  costs  for  the  three  approaches 
(and  it  must  be  realized  that  because  no  design  detail  is  available,  these 
estimates  are  very  rough): 


Software 

Hardware 

Design 

Systems 

Total 

Equivalent 

Cards 

GPSP 

200  m-day s 

120  m-days 

300  m-days 

620  m-days 

80 

UFFT/CPE 

200 

120 

300 

620 

80 

CE 

200 

500 

400 

1100 

135 

We  can  combine  the  variable  hardware  costs  with  the  fixed 
development  costs  to  form  the  curves  of  Figure  5-13  and  5-14.  Note  that 
the  CiPSP  looks  more  attractive  here,  because  of  the  higher  development 
cost  associated  with  the  CE  approach. 


5-33 


1 


K**S**>  xWy  + 


Comparison:  Hardware  Plus  Development  Costs 


igure  5-14.  Performance  Available  for  Constant  Total  Cost  (C  = 200) 


SECTION  6.  CONCLUSIONS  AND  RECOMMENDATIONS 


6.  1 The  elements  of  a CFAR  signal  processor  have  been  defined  and 
described  in  some  detail.  It  is  clear  that  the  coherent  processing  offered 
by  the  FFT  will  improve  radar  sensitivity,  and  that  the  CFAR  circuitry 
suggested  here--range-cell  averaging,  clutter  mapping- -should  suppress 
common  types  of  false  target  sources. 

6.2  The  limits  on  processing  will  be  set  ultimately  by  cost.  Some 
insights  have  been  gained  into  the  relative  costs  of  various  processing 
approaches,  as  well  as  the  relationship  between  cost  and  performance 
parameters  such  as  number  of  pulses  coherently  integrated  and  number  of 
range  cells  handled. 

6.3  Further  refinement  of  the  analysis  would  be  desirable.  This  is 
particularly  true  in  the  following  areas: 

(a)  MTI:  We  have  investigated  only  two  MTI  approaches:  three-pulse 

MTI  and  narrow-notch  N-pulse  MTI.  Neither  is  felt  to  be  optimum  for  the 
EAR  environment. 

(b)  Clutter-map  parameters:  More  data  would  be  desirable  on  the 

effect  of  clutter-map  parameters  (i.  e. , the  number  and  "age"  of  clutter 
samples  averaged  and  stored)  on  system  performance. 

(c)  Quantization:  Although  the  present  results  are  good  enough  to  be 

used  with  the  present  9-bit  EAR  A/D  converter,  no  analysis  exists  to  tell 
us  whether  anything  would  be  gained  if  (at  some  future  time)  it  were  decided 
to  increase  the  number  of  bits  in  the  processor. 

(d)  Coherent  versus  Non-Coherent  Integration:  For  a given  beam 

position  dwell  time,  there  is  an  interesting  trade-off  to  be  examined  between 
coherent  integration  (FFT)  and  incoherent  integration,  with  the  possibility 

of  rf  frequency  changing  between  coherent  pulse  groups  taken  into  account 
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as  a means  of  reducing  target  fluctuation  losses.  This  tradeoff  affects 
sensitivity;  the  effect  of  such  integration  changes  on  CFAR  performance, 
which  is  more  complex,  also  needs  to  be  taken  into  account. 

6.4  In  the  hardware  area,  it  is  clear  that,  although  a substantial  beginning 
has  been  made,  the  level  of  detail  is  insufficient  for  a design  to  be  initiated. 
But  we  have  identified  three  usable  processor  techniques  (CE,  GPSP,  and 
UFFT).  The  Common  Element  (CE)  approach  is  the  newest,  and  as  a result 
entails  more  development  cost.  B ut  it  has  the  promise  of  least  ultimate 
hardware  cost,  as  well  as  being  an  inherently  modular  design,  with  great 
flexibility  in  terms  of  both  re-programmability  and  growth  to  more  range 
cells/samples.  A development  program  can  be  visualized  which  starts  out 
modestly  with  the  development  of  only  one  Common  Element  module,  upon 
which  various  processing  algorithms  can  be  exercised.  Ultimately,  the 
modules  can  be  replicated  in  sufficient  quantity  to  form  a total  processor  of 
whatever  capability  is  desired. 
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PROGRAMS  AND  FLOW  CHARTS 
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Program  CFAR  Page  1 of  12 


t 

r 

i 


I 


V 


PROGRAM  CFAR  < INPUT  t OUTPUT  » TAPES- INPUT » T APE6-0UTPUT » TAPE2 > 

C 

C PROGRAM  CFAR  WRITTEN  FOR  P. CORNWELL  GY  C. 8. COREY,  APRIL  1976 
C 

c 

C INPUTS 
C 

C NPT8-NUMBER  OF  POINTS  IN  FFT 
C RHOG-GROUND  CLUTTER  CORRELATION  COEFFICIENT 
C NRCELLS-NUMBER  OF  RANGE  CELLS 
C UMDL -WEATHER  MODEL  *1  OR  WEATHER  MODEL  *2 
C NT I -TURNS  MTI  FILTER  ON  OR  OFF 
C MTIMDL-MTI  MODEL  *1  OR  MTI  MODELA2 
C THMDL-TNRE8HOLD  MODEL  *1  OR  THRESHOLD  MODEL  92 
C NCLMAP-NUMBER  OF  POINTS  IN  CLUTTER  MAP 
C 8TNDB-8IGNAL  TO  NOISE  RATIO  IN  DG 
C STGCDB-8 I GNAL  TO  GROUND  CLUTTER  RATIO  IN  DG 
C A- WEI BULL  SLOPE  PARAMETER 
C T 

C IWGT-TURN8  WEIGHTS  ON  OR  OFF 
C S I G8G-3PECTR AL  WIDTH  OF  GROUND  CLUTTER 
C 8I08W-8PECTRAL  WIDTH  OF  WEATHER  CLUTTER 
C GH-TURNS  GROUND  CLUTTER  MODEL  ON  OR  OFF 
C NM-TURN8  NOISE  MODEL  ON  OR  OFF 
C WM-TURN8  WEATHER  MODEL  ON  OR  OFF 

C MUW1-MEAN  WEATHER  SIGNAL  TO  LEFT  OF  CELL  OF  INTEREST 
C MUW2-MEAN  WEATHER  SIGNAL  TO  RIGHT  OF  CELL  OF  INTEREST 
C FW-CENTER  FRQUENCY  FOR  RAIN  STORM 
C MUWO-MEAN  WEATHER  SIGNAL  CELL  OF  INTERE8T 
C KCLUT-C0N8T ANT  CLUTTER  FLAG 
C 8TKCLDB-SIGNAL  TO  CON8TANT  CLUTTER  IN  DG 
C NBAR-PARAMETER  FOR  TAYLOR  WEIGHTS  CALCULATION 
C ATAY-PARAMETER  FOR  TAYLOR  WEIGHTS  CALCULATION 
C ITATTN-PLOT  OPTION  FOR  CLUTTER  MAP  THRESHOLD  ATTENUATION  V8. 

C FREQUENCY  PLOT 

C IATTN-PLOT  OPTION  FOR  GROUND  CLUTTER  ATTENUATION  MS.  FREQUENCY  PLOT 
C MORCFAR- 
C 
c 

LOGICAL  MTI'QM'WM'KCLUT'IWGT'NMrITATTN'IATTN'MORCFAR 
C 

REAL  MUT , NUN  t MUCL  * KCLT8 » KCL  r HUCLF  t MUWO  » HUN 1 » MUW2 
C 

INTEGER  CHECK  t WMDL » THMDL  * WFL AG 
C 

DIMENSION  WGTS(64>'R<66'66)'B(66'66>'PZ(66'66>' 

1 GT<66'66> >PZTEMP<66'66) »PZTIL<64r64>  »C(64»64> »CT<64r64> » 

2 PW(64f 64) rCP(64,64> ,D(64,64) »8IGC8Q(64) »8IGW8Q(64) »8IGNSQ<64> 


ii  i ■ ■ ^ -- «*■,«  % - --  - - 
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c 

EQUIVALENCE  (CTfBTfCPfATTNfTATTNfTI fSTIDB) f (DfRf ATTNDBfTATTNDB) 
EQUIVALENCE  <BfPU> 

C 

DATA  PI/3.14150245/»T/2.0E-04/ 

DATA  RHOG  » NRCELL8 » MMDL t NT I / 1 ♦ 0 1 4 1 1 » . T . / » NT I MDL r THNDL » 

1 NCLMAPrSTNDB/1 t 1 f4» 10.0/f8TGCDBf Af IWGT f8IG8G/-35»0f 

2 2.0f .T. f 18.0/f8IGBWfGMfNMfWN/108*0.f  .T.r.T.r.T./r 

3 HUW1 »HUW2fFWfMUW0/3162»0f3142.0f25O0*0f3162.0/»KCLUT/.F./» 

4 ITATTNf IATTNt  MORCFARAF*  t *F» t*F*/ 

C 

NANEL I ST / 1 NP3/HUW 1 » MUW2 » MUWO f FW 
NAHEL I 8T / I NP2/8TKCLDB 

NANEL I ST/ INP 1 /NPTS  » WNDL r NT I » NT IMDL  r THNDL  » 

1 STNDBfSTOCDBf A» IM6T»81Q86r8IQ8N» ©NrNNfNNfKCLUT fNBARfATAY 

C 

C INPUTS 
C 

PRINT  801 

801  F0RMAT</f10Xf*INPUT  XNP1  - NPTSfWNDLfMTI f*f// , 

1 10Xf*MTIMDLfTHMDLfSTMDBfSTGCDBfAfIWGTfSIGSGf*f//f 

2 10XF*8IQ8Wf 8NfNN»UHrKCLUT»NBAR»ATAYBf /> 

READ  INP1 

PRINT  INP1 
C 

C INITIALIZATION 

C 

CHECK-0 
UFLAG-0 
DO  6 I-lfNPT8 
SI0C8Q( I )-0. 

8I0W8Q( X )-0* 

SI0N8Q(I)-0. 

WGT8( I )-l * 

6 CONTINUE 
C 

C CONPUTE  NEAN  OF  TARGET f NOISE  AND  CLUTTER 
C 

NUT-1 .0 

NUN-NUTB 1 0 . OB* < -8TNDD/ 1 0 . 0 ) 

NUCL-NUTB10 . OB* ( -8T8CDB/10 . 0 > 

C 

C 18  NT!  MODEL  TO  BE  CONPUTEDT 
C 

IF < *NOT *NTX  >00  TO  50 


f i ; c 

[f  C IF  YES t WHICH  MODEL  18  TO  BE  U8EDf  1 OR  2T 


Program  CFAU  (Continued) 
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NPTSP2-NPTS+2 
8108080-81080**2 
T80-T**2/2»0 
DO  19.  X-1»NPT8P2 
DO  13  J-lrNPT8P2 

R< I» J)-EXP( <-<2.0*PI*8I080*(I-J)*T)**2)/2«0> 

15  CONTINUE 

C********** 

C OUTPUT  CR<X*J>3  FOR  OROUND  CLUTTER  (NPT8+2  X NPTS+2) 

C********** 

C PRINT  P08 

908  FORMAT <//* IX *»  R MATRIX  FOR  OROUND  CLUTTER**/ > 

DO  14  1-1»NPT8P2 

C PRINT  902*  <R< I » J) * J-l * NPT8P2) 

14  CONTINUE 
C 

C INITIALIZE  CD< X * J)3  MATRIX  FOR  MODEL  1 (NPT8+2  X NPT8+2) 

C 

18  NPT3P2-NPT8+2 
DO  20  I-1*NPT8P2 

DO  20  J-l »NPT8P2 

B( X * J)— 0*0 
20  CONTINUE 

DO  23  1-1 rNPTB 

b<x»x)-i;o 

B(I*I+1>— 2.0 
B<  I » I+2)-l *0 
25  CONTINUE 

B ( NPT8P2- 1 * NPT8P2- 1 ) - 1 . 0 
B ( NPT8P2  » NPT8P2 ) - 1 .0 

C 

C COMPUTE  CB<I* J)3  X CR<X*J>3  <NPT8+2  X NPT8+2) 

C 

DO  30  1-1 *NPT8P2 
DO  30  J-l *NPT8P2 
PZTEMP< I * J>-0.0 
DO  30  K-l » NPT8P2 

PZTEMPf  X* J)-PZTEMP<  X* J>+B(X*K>*R(K* J> 

30  CONTINUE 
C 

C COMPUTE  CBT(I* J)3 
C 

DO  35  I-i»NPTSP2 
DO  33  J-1*NPT8P2 
BTd  » J)-B(  J*I ) 

33  CONTINUE 
C 

C COMPUTE  CPZd*  J)3-CBd*  J)3*CRd*  J>3*(BTd*  J)3  (NPT8+2  X NPT8+2) 

C 

DO  40  I-1*NPT8P2 
DO  40  J-1*NPT8P2 
PZd*J>-0.0 
DO  40  K-1*NPT8P2 

PZ< I » J)-PZ< I » J)+PZTEMP( X *K)*BT  <K*  J) 

40  CONTINUE 
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C 

C OUTPUT  CPZd*J>3 

c 

C PRINT  901 

901  FORMAT <//» IX *«  PZ  MATRIX  HTZ  MODEL  !*»/) 

DO  42  X-1»NPTSP2 

C PRINT  902»<PZ(X'J)'J-1»NPT8P2> 

902  F0RMAT</»1X'4E13.4»/»10UX'4E13.4'/)) 

42  CONTINUE 

C 

00  TO  40 
C 

C NO  MTI 
C 

C 18  CONSTANT  CLUTTER' TO  BE  COMPUTED? 

C 

SO  IF ( .NOT »KCLUT)GO  TO  S7 
C 

C YES-INPUT  SIGNAL  TO  C0N8TANT  CLUTTER  IN  DB  AND  COMPUTE  C0N8TANT  CLUTTER 
C 

PRINT  802 

802  F0RMAT(/'10X'«XNPUT  XNP2  - STKCLDBt » / ) 

READ  INP2 
PRINT  INP2 

KCL  T8» 1 0 » 0* * ( -8TKCLDB/ 1 0 . 0 ) 

KCL-MUTBKCLT8 

MUCLF- 1 0 . 0**  < -8T0CDB/ 10.0) -KCLT8 

C 

C TEMPORARY  STOP  - A* 

C 

CALL  EXIT 
C 

C NO  - COMPUTE  CRd»  J)3  AND  CBd*J)3  MATRXCE8  CNPT8  X NPT8) 

C 

37  SIGSGSQ-8IG8G**2 
TSQ-TM2/2.0 
DO  43  I-1.NPT8 
DO  43  J-1'NPTB 

Rd'J)«EXP<<-<2.0*PX*SXGSG*d-J)*T)**2)/2.0) 

Bd*J)>0.0 

43  CONTINUE 

DO  47  I-l.NPTS 
DO  47  J-lrNPTS 
IFd.NE. J)00  TO  47 
Bd»J)-l*0 
47  CONTINUE 
C*«««»««**« 

C OUTPUT  CRd'J>3  NO  C0N8TANT  CLUTTER 
C*««*B****« 

C PRINT  909 

909  FORMAT <//» IX »«  NO  CONSTANT  CLUTTER  R MATRIX**/) 

DO  44  X-lrNPTS 

C PRINT  902'<Rd»J)'J>l*NPT8) 

44  CONTINUE 
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C 

CHECK-1 
00  TO  AO 

C 

C HTI  MODEL  2 
C 

C COMPUTE  INITIAL  COVARIANCE  MATRIX  FOR  GROUND  CLUTTER  (NPT8  X NPT8) 
C 

43  8168080-81080**2 
T8Q-T«*2/2 « 0 
DO  69  I-ltNPTS 
DO  69  J-l fNPTS 

R( Ip J)-EXP( (-(2.0*PI*8I0B0*( I-J)*T)**2)/2»0) 

69  CONTINUE 

C********** 

C OUTPUT  CR(Ip J)3  MTI  MODEL  2 
C********** 

C PRINT  910 

910  FORMAT (//p IX p*  R MATRIX  MTI  MODEL  2*p/> 

DO  71  1-1 fNPTS 

C PRINT  902p(R(IpJ)pJ-1pNPT8) 

71  CONTINUE 
C 

C DEFINE  CB( I p J) 3 FOR  MODEL  2 (NPT8  X NPT8) 

C 

130  DO  74  1-1 fNPTS 
DO  74  J-l pNPTS 

IF( I »EQ» J)B( I r J>— ( 1 *0-1 •O/NPTB) 

IF(I»NE»J)B( I p J)— 1 *0/NPT8 
74  CONTINUE 
C 

CHECK- 1 
C 

60  IF < CHECK «NE* 1)00  TO  70 
C 

C COMPUTE  CPZTEMP(IpJ)3-CB(IpJ>3*CR(IpJ>3  FOR  MTI  MODEL  2 OR  NO  MTI 
C (NPT8  X NPTS) 

C 

DO  75  1-1 fNPTS 
DO  73  J-1fNPT8 
PZTEMP(IpJ)-0»0 
DO  75  K-1pNPT8 

PZTEMP  < I f J ) -PZTEMP (IfJ)+B(I»K)*R(K»J) 

73  CONTINUE 

C 

C COMPUTE  CBT(IfJ)3  FOR  MTI  MODEL  2 OR  NO  MTI  (NPTS  X NPTS) 

C 

DO  80  1*1 f NPTS 

DO  80  J-l fNPTS 
BT(XfJ)— B(JfI> 

80  CONTINUE 

C 

C COMPUTE  COVARIANCE  MATRIX  CPZTILUp J)3-CB(If J)3*CR(If J)3*CBT< If J>3 
C FOR  MTI  MODEL  2 OR  NO  MTI  (NPTS  X NPTS) 
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c 

DO  85  I-1»NPT6 

DO  85  J-1'NPTS 

PZTIL<I»J>-0.0 
DO  85  K-1'NPT8 

PZTILflt J>»P2TIL<I* J)+PZTlHP<X»K)8tT(K» J) 

85  CONTINUE 
Cttttttttlt 

C OUTPUT  CPZTILU»J>3  FOR  NTX  MODEL  2 OR  NO  MTI 
C***tSt**** 

C PRINT  903 

903  FQRMAT<//»lXr*  PZTXL  MATRIX  MTI  MODEL  2 OR  NO  MTX8»/> 
DO  91  I«1»NPT8 

C PRINT  902 » (PZTILdt  J>  * J*1»NPT8> 

91  CONTINUE 
C 

00  TO  90 
C 

C COVARIANCE  MATRIX  AFTER  MTI  MODEL  2 CPZTIL<I»J>3  (NPT8  X NPT8) 
C 

70  NPT8P2-NPT892 
DO  72  X-1'NPTS 

DO  72  J«1»NPT8 

PZTXL<If J)-PZ<IfJ) 

72  CONTINUE 
Ct«****t**t 

C OUTPUT  CPZTILd»J)3  AFTER  MTI 
CS*8««*S*** 

C PRINT  913 

913  FORMAT <//» IX »*  PZTXL  MATRIX  AFTER  MTX*»/> 

DO  93  I»1#NPTS 

C PRINT  902f <PZTXL< I f J) » J*lt NPT8) 

93  CONTINUE 
C 

C COMPUTE  HEIGHTS  IF  DESIRED 
C 

90  IF < »NOT • I WOT ) GO  TO  95 
IF(WFLAG.EQ.1)G0  TO  98 
CALL  TAYLOR ( HOTS fNPT8»NBARrATAY> 


WFLAG-1 
Cf  ******** 

C OUTPUT  HEIGHTS 
CS*»*****S» 

98  CONTINUE 
C PRINT  921 

921  FORMAT <//» IX »t  TAYLOR  HEXQHT8t'/> 

DO  149  I«1 »NPT8 

C PRINT  922rXrHGT8<X) 

922  FORMATdXr X5fE15*8> 

149  CONTINUE 


C 

C DEFINE  HEIGHT  MATRIX  CC<X'J>3 


<NPT8  X NPT8> 


1 


1 


f 
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c 

DO  92  I«1*NPTS 
DO  92  J-1*NPT3 
C<I*J>-0.0 

IF<  I »EQ» J)C(I * J)"WGT8( I > 

92  CONTINUE 

C 

C DEFINE  COVARIANCE  MATRIX  AT  FFT  INPUT  CPWCI* J)>CC<I* J>1*CPZTIL(I* J>3 
C «CCT<I*J>3  (NPT8  X NPT8) 

C 

DO  94  1*1 *NPTS 
DO  94  J-1'NPTS 
CT(I*J)-C<J*I> 

PZTEMPC I * J)*0»0 
DO  94  K-1'NPT8 

PZTEMPtI* J)-PZTEMP<I* J)+C< I *K>*PZTIL<Kr J) 

94  CONTINUE 

DO  96  I«1 *NPTS 
DO  96  J-lrNPTS 
PU(IrJ)*0.0 
DO  96  K“1 t NPT8 

PW ( I * J ) «PW ( I * J ) +PZTEHP  < I * K > *CT  ( K * J ) 

96  CONTINUE 
C********** 

C OUTPUT  CPM<I* J>3 
C********** 

C PRINT  920 

920  FORMAT*//* 1X»*  PM  MATRIX  WITH  MEIGHT8  AT  FFT  INPUT**/) 

DO  97  I-1*NPT8 

C PRINT  ?02»<PM<I'J)'J-1'NPTS> 

97  CONTINUE 
C 

00  TO  100 
C 

C NO  WEIGHTS 
C 

93  DO  99  I-1*NPT8 
DO  99  J“1 *NPTS 
PU<I*J>-PZTIL<I» J) 

99  CONTINUE 
C 

100  CONTINUE 

IF  < I WOT >00  TO  102 
C*****t*tt* 

C OUTPUT  PW 

c*******t«t 

C PRINT  904 

904  FORMAT <//* IX »*  PW  MATRIX  WITHOUT  WEIGHTS**/) 

DO  101  I«1*NPT8 

C PRINT  902* <PW(I* J) * J"1*NPT8) 

101  CONTINUE 


C 

C 

C 


IS  GROUND  CLUTTER  TO  BE  COMPUTED? 


102  IF<  *MOT .GM>GO  TO  103 
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C YES  - COMPUTE  CCP(X*K)3  AND  CD(X*K)3  FOR  0ROUND  CLUTTER  MODEL 
C 

DO  110  I-ltNPTS 
DO  110  K-1*NPT8 
KU8E-K-1 

CP(I*K)-C0S<<2.0*PX)/NPT8*(I~1>*KU8E) 

D< I*K)-8IN< (2.08PX )/NPT8*( X-l )*KU8E) 

110  CONTINUE 
C********** 

C OUTPUT  CCP(I* J)3  AND  CD(I*J)3  FOR  GROUND  CLUTTER  MODEL 
C********** 

C PRINT  914 

914  FORMAT <//* IX**  CP(I» J)  MATRIX  FOR  GROUND  CLUTTER** /> 

DO  141  1-1 *NPT8 

C PRINT  902*  <CP<I * J) * J-1*NPTS) 

141  CONTINUE 

C PRINT  915 

915  FORMAT*//* IX »*  D(IrJ)  MATRIX  FOR  GROUND  CLUTTER**/) 

DO  142  I-1*NPTS 

C PRINT  902*<D(I* J) * J-1*NPT8> 

142  CONTINUE 

C 

C COMPUTE  DOPPLER-DEPENDENT  GROUND  CLUTTER  VARIANCE  MODULO  RANDOM 
C CONSTANT  WHICH  IS  RANGE  DEPENDENT  (EACH  CHANNEL) 

C 

DO  115  K-l *NPT8 
SUM-0. 0 

DO  116  I-l *NPT8 
DO  116  J-l *NPT8 

8UM-SUM+PW< I * J)*CP( 1 »K>*CP< J»K>+PW< I *J)*D<I*K)*D(J*K) 

116  CONTINUE 

SIGCSOOO-SUM 
115  CONTINUE 
C********** 

C OUTPUT  SIGCSQ 
C********** 

PRINT  905 

905  FORMAT <//> IX **  8IGC8Q  ARRAY**/) 

PRINT  902* (8IGC8Q(I ) *X>1*NPT8) 

C 

C SET  GROUND  CLUTTER  MODEL  FLA0-FAL8E 
C 

OM-.F. 

C 

C 18  WEATHER  CLUTTER  MODEL  TO  BE  COMPUTED? 

C 

IF(WM)GO  TO  125 
C 

C NO  - REINITIALIZE  WEATHER  PARAMETERS 
C 

8IG8W-0.0 

MUW1-0.0 

MUW2-0.0 

MUWO-O.O 

FW-0.0 

DO  103  I-l *NPT8 
8IGW8Q( I )-0.0 
103  CONTINUE 


i 


i 
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C NO  WEATHER  TO  BE  COMPUTED  - 00  TO  N0X8E  MODEL 

C 

00  TO  120 
C 

C GROUND  MODEL  18  NOT  TO  BE  CALCULATED 
C 

105  IF( .NOT »WM)0O  TO  130 
C 

C COMPUTE  CCP<I»K)3  AND  CD<X'K)3  FOR  WEATHER  CLUTTER  MODEL 
C 

DO  121  I-lrNPTS 
DO  121  K-1'NPTS 
KUSE-K-1 

CP< I *K )«COB< 2 »0*PI* < 1-1 ) <<FL0AT<KU8E )/NPTS)  > 
D<IrK)-8IN<2.0*PI*<I-l)*<FL0AT<KUSE)/NPTS)> 

121  CONTINUE 
C*******«*t~ 

C OUTPUT  CCP< X r J>3  AND  CD<1'J>3  FOR  WEATHER  CLUTTER  MODEL 
C********** 

C PRINT  916 

914  FORMAT < //r IX r*  CP<IrJ)  MATRIX  FOR  WEATHER  CLUTTER**/) 

DO  143  X-l »NPTB 

C PRINT  902*  <CP<I» J)  * J-1*NPTB> 

143  CONTINUE 
C PRINT  917 

917  FORMAT <//* IX **  D<X*J>  MATRIX  FOR  WEATHER  CLUTTER* * / ) 

DO  146  1-1 *NPT8 

C PRINT  902*<D<I*J)* J-1*NPT8> 

146  CONTINUE 

C 

C COMPUTE  DOPPLER  DEPENDENT  WEATHER  CLUTTER  VARIANCE  MODULO  A CONSTANT 
C WHICH  18  RANGE  DEPENDENT  (EACH  ELEMENT) 

C 

DO  122  K-1'NPTS 
SUM-0.0 

DO  124  1-1 *NPTS 
DO  124  J— 1 »NPTB 

8UM-8UM+PW< I * J)*CP< I *K)*CP< J*K)+PW< I * J)*D< I *K)*D< J*K) 

124  CONTINUE 

SX0W80(K)-8UM 

122  CONTINUE 
C********** 

C OUTPUT  8X8W80 
C*******»ft* 

PRINT  906 

906  FORMAT!//* 1X*«  8X0W8Q  ARRAY**/) 

PRXNT  902*  <8IGW80< I ) * X-l *NPTB) 

C 

C SET  WEATHER  CLUTTER  MODEL  FLAG-FALSE 
C 

WM-.F. 

00  TO  120 
C 

C NO  WEATHER  CLUTTER  - COMPUTE  NOISE  CLUTTER 

C 
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C COMPUTE  CCP(1*K)3  AND  CD(X*K)3  PON  NOISE  CLUTTER 
C 

130  DO  131  X«1«NPT8 
DO  131  K«1 *NPT8 
KU8E-K-1 

CP < I * K > -COS < < 2 . 0*PI ) /NPT8* < X-l ) SKU8E ) 
D<I*K)"8XM<<2.0*PX)/MPT8*<X-1)*KU8E) 

131  CONTINUE  ' 

Cttttttlttt 

C OUTPUT  CCP(I*J)3  AND  £D(X*J)3  FOR  N0X8E  CLUTTER  MODEL 
C********** 

C PRINT  918 

918  FORMAT (//* IX **  CP(X*J)  MATRIX  FOR  NOISE  MODEL**/) 

DO  147  I-1*NPTS 

C PRINT  902*  <CP(X* J) * J*1»NPT8> 

147  CONTINUE 

C PRINT  919 

919  FORMAT <//, IX**  D<I*J>  MATRIX  FOR  N0I8E  MODEL**/) 

DO  148  I-1*NPT8 

C PRINT  902*<D(X*J)*J-1*NPT8) 

148  CONTINUE 
C 

C COMPUTE  DOPPLER  DEPENDENT  N0I8E  CLUTTER  VARIANCE  (EACH  CHANNEL) 
C 

DO  132  K-1*NPTS 
8UM-0.0  * 

DO  133  I“1 *NPT8 
DO  133  J-1*NPT8 

8UM"8UM+PW(  I * J)«CP(  X *K)*CP( J*K)+PW(  I*J)*D(I*K)«D(J*K) 

133  CONTINUE 

8IGN80(K)-8UM 

132  CONTINUE 
C********** 

C OUTPUT  8I0N80 
C********** 

‘ ‘ ""PRINT  907 

907  F0RMAT(//*1X*«  8X8N8Q  ARRAY**/) 

PRINT  902* (8X8N8Q(X) *I"1*NPT8> 

MR I TE < 6 * 1 000 ) HUT * MUCL * KCL * NUN * MUCLF * 

+MUW1 * MUU2 » MUWO 

1 000  FORMAT  < 6X*MUT-*E20 . 8//6X*MUCL-*E20  * 8//6X* 

+*KCL-*E20 . 8//6X*MUM-*E20 . 8//6X  * 

♦*MUCLF-*E20 . 8//6X  * *MUM 1 -*E20 . 8//&X * 

+*MUW2-*E20 . 8//6X * *MUWO-*E20 , 8 ) 

MR I TE ( 2 ) HUT » MUCL * KCL * MUN * MUCLF * MUNI  * MUN2 * MUNO 
HRITE(2)(8I0C80<K)*K-1*NPT8)*(8I0M8Q(K)*K-1*NPT8)* 
+<8I0N80<K) *K-1*NPT8> 

MR1TE(2) CMOTSOO  *K«1 rNPTS) 

ENDFILE2 

C 

00  TO  140 
C 

C WEATHER  CLUTTER  MODEL 
C 

C HEATHER  MODEL  1 OR  MODEL  2T 
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C 

125  IF(WMDL.E0»2>CALL  EXIT 
C 

C WEATHER  MODEL  1 
C 

PRINT  003 

803  FORMAT </r*INPUT  IMPS  - MUWlfMUW2rMUW0rFW«> 

READ  I HP 3 
PRINT  INP3 

C 

C COMPUTE  CR( I » J) 3 MATRIX  FOR  WEATHER  MODEL  1 
C 

NPT8P2-NPT8+2 
T8Q-T**2/2.0 
8IB8W8Q-8IG8Wt*2 
DO  135  I»1 r NPTSP2 
DO  135  J-1»NPT8P2 

RUtJ>-EXP<<-<2.0*PI*8IG8W*d-J>*T>**2>/2.0> 

Rdf  J>-RdrJ>8C08<2.0«PI*FW*d-J>*T> 

135  CONTINUE 
Cttf******* 

C OUTPUT  CRdf  J)3  FOR  WEATHER  MODEL  1 
C******ttft 
C PRINT  912 

912  FORMAT <//, IX »*  R MATRIX  WEATHER  MODEL  1 *»/> 
DO  134  I»1 f NPTSP2 

C PRINT  902r  <R< I » J) » J"1 »NPT8P2> 

134  CONTINUE 

C 

C CHECK  IF  HTI  BEING  DONE 
C 

144  IFtMTDGO  TO  145 
C 

C NO  HTI  - DEFINE  CBdrJ>3 
C 

DO  138  I-1'NPTS 
DO  138  J-1»NPT8 
B(Ir J)-0.0 
IFd.EQ. J)B<Z»J>-1.0 
138  CONTINUE 
C 

GO  TO  60 
C 

C CHOOSE  MTI  MODEL  1 OR  MODEL  2 

C 

145  IF<MTIMDL.E0.2>G0  TO  150 
GO  TO  18 

C 

C NO  WEATHER  CLUTTER  OR  WEATHER  CLUTTER  ALREADY  DEFINED 

C 

C N0I8E  MODEL 
C 

120  NPT8P2-MPTB+2 

DO  123  I-1'NPT8P2 
DO  123  J-lfNPT8P2 
R(lf J>-0.0 
IFd.E0.J>R<IfJ)«l»0 
123  CONTINUE 
GO  TO  144 
C 


1 

}i 

i 
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C COMPUTE  CLUTTER  MRP  THRESHOLD  ATTENUATION  MiJ  CONVERT  TO  DB 
C 

140  CONTINUE 
CALL  EXIT 
END 

C«««****t*«**«*««»*»*B««*t**«m»««**t«*M»*«***tS««t****t*tt»tt** 
SUBROUTINE  TAYLOR<GT»NGT*NBAR»A> 

C 

C COMPUTATION  OF  THE  TAYLOR  ILLUMINATION  FUNCTION 
C GT-ARRAY  TO  BE  RETURNED  TO  MAIN  PROGRAM  CONTAINING  TAYLOR  HEIGHTS 
C NGT-DIMEN8I0N  OF  GT 
C 

DIMENSION  GT<NGT)»F<64> 

PI-3.14139243 
NBARM1-NBAR-1 
FMBAR— NBAR 

SIG80-FNBAR*»2/ ( A«*2+  < FNBAR- . S ) ««2 ) 

C 

N0T2-N0T/2 

FNGT-NOT 

DO  100  IXPR-1»NGT2 
XPR-IXPR 

XPRlME-< ( (-FNGT/2.0) )+XPR)/FNGT 

8UMFM 1-0.0 

SUHFM2-0.0 

DO  SO  M-lrNBARMl 

FM-M 

PROD 1-1.0 

PR0D2-1.0 

DO  40  N-lrNBARMl 

FN-N 

ARG1 A-A**2+ ( FN- . 3) **2 
ARG1-1 .0-FM**2/(8IG8Q*AR01A> 

PR0D1-PR0D1BARG1 
IF<N«EQ.M)GO  TO  40 
ARG2- 1 . 0- < FM**2/FN**2 > 

PR0D2-PR0D2BARG2 
40  CONTINUE 

F<M)-<(<-l . 0)**(N+1 ) ) SPR0D1 )/(2.0tPR0D2) 

SUMFH 1 -8UMFM 1 +F  < M ) *C08 ( 2 . OtPI BFMtXPRI ME  > 

8UMFM2-SUMFN2+F  < M > 

30  CONTINUE 

GT < IXPR > - < 1 . 0+2 . OtSUMFHl ) / < 1 . 0+2 . 0S8UHFH2 > 

100  CONTINUE 
C 

DO  200  I-1»NGT2 
IM1-I-1 

GT(N0T2+I >-GT<NGT2-IMl > 

200  CONTINUE 

C 

RETURN 

END 

ctiiiTitiiivfmmmmmmtmmtmmmmtmmmm 

READY » 


R333114  LOG  OFF 
R333114  8RU 


10.49.03. 
2.363  UNTS. 
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PROGRAM  CFARP< INPUT » OUTPUT » TAPES- INPUT ► TAPE4-0UTPUT » TAPE 2 ) 
REAL  MUT « MUCL » KCL  » NUN  » MUCLF » MUNI » HUU2 1 HUWO 
INTEGER  THMDL' RHOG 

LOG I CAL  I TATTN  f I ATTN  f XFRE3P  r MORCFAR » PRF A » I DEALCF 
LOGICAL  STIPL»MTI»IUGT tPRED 

DIMENSION  8IGCSQ<64)  »SIGM8Q<64)  »8IGN8Q<64> »TI<64) 

DIMENSION  Til (64) »8TX1<64>  »CMT<64>  »CFAR1<64) 

DIMENSION  CFAR2(64> » BETA < 44) tPFA(64> *PFA6UM<64> 

DIMENSION  TTE8T<64) »E<64) 'PD<44) »8UMPD(64) 

DIMENSION  STI<64) '8TXDB(64>  »8CALET<64) »WGTS<64>  »SIG< 1000) 
DIMENSION  T AT TN  < 64  > » TATTNDB ( 44  > » FREQ ( 64  > r ATTN < 64  > r ATTNDB ( 64  > 
DIMENSION  FRE8P(64)rFRE8PDB<64> 

DATA  T'PI/2.E-04t 3. 14159265/ 

DATA  A/2./ 

DATA  NUU/1000/ 

NAMEL I 8T/0UT 1 /S I G 

NAMEL 1 8T/0UT2/CF AR 1 1 CFAR2 » BETA  t PFA  t PFA8UM » TTE8T  »E.PD» 8UMPD t 
1CMT 

NAMEL IST/0UT3/HUT  » MUCL  » KCL • MUNr MUCLF » MUNI r MUW2.MUW0 
NAMEL 1 8T/0UT4/8  XGC80 1 8 X QU8Q r 8XGN8Q t WOTS 
NAMEL I 8T/0UTS/8C ALET 

NAMELIST/ INP1 /RHOG  t NRCELL8 1 NCLHAP t NPTS » MTXMDL 

1 ITATTN  » I ATTN  » IFRE8P » MORCFAR  t MTI » IUGT  r 8TIPL r THMDL • 

2 PRFAf I DEALCF rPRED 
WRITE(AflO) 

10  FORMAT <6X*INPUT  I NP 1 1 RHOG » NRCELL8 1 NCLMAP r NPTS  t MT IMDL  » * r / 

1 tlTATTNt I ATTN . XFRE8P  * MORCFAR  » MTI » IUGT  » 8TIPL  r THMDL*/ 

2 *PRFA» I DEALCF »PRED*) 

READ  XNP1 

REUIND2 

READ < 2 ) MUT » MUCL r KCL » MUN r MUCLF t HUU1 » MUU2 r MONO 
READ < 2 X 8XGC8Q < K ) » K-l t NPTS > r < 8X0W8Q ( K ) » K-l » NPTS > » 

1 (SIGNSQ(K) »K-1 » NPTS) 

READ(2) (WGT8(K> 'K-l »NPT8> 

140  DO  153  K-l 'NPTS 

IF (SIGC8Q<K) .LT.O* >SIGCSQ(K)-1 .E-50 
IF  <8IGU8Q(K> .LT.O. )SXGU8Q(K>-1 .E-50 
IF(SIGN8Q<K) .LT.O. >SX0N8Q(K>-1 .E-50 
TEST-MUCL*SIGCSQ( 1 > 

IF < TEST . EG . 0 . 0 > TEST- 1 . E-50 
TATTN<K)-(MUCLtSXGC8Q(K>  >/ 

1TEST 

TATTNDB  < K ) - 1 0 . 0* AL0G1 0 ( TATTN  <K>  > 

IF  < TATTNDB(K) .LT. (-100.0) > TATTNDB <K>— 100.0 
153  CONTINUE 

C 

C CHECK  IF  CLUTTER  MAP  THRESHOLD  ATTENUATION  VS. FREQUENCY  IS  TO 
C BE  PLOTTED 
C 

IF< .NOT. I TATTN) GO  TO  155 

C 

C COMPUTE  FREQUENCY  ARRAY 
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FACT-1. 0/<T*HPTS> 

DO  140  K-1.NPT8 
FREQ<K)-<K-1)SFACT 
140  CONTINUE 

~ C 

C PLOT  CLUTTER  HAP  THRESHOLD  ATTENUATION  V8. FREQUENCY 
C 

PRINT  923 

923  FORMAT < // » 1 X » SCLUTTER  HAP  THRESHOLD  ATTENUATION  VS. FREQ*  PLOTS'/* 

1 lX'SSET  UP  PAPER  FOR  0 TO  -BO  DB*  AND  TYPE  60S) 

READ  804 'IQO 
804  F0RMAT<A2) 

CALL  PLOT < FREQ »TATTNDB»NPTS> 

C 

“ C PRINT  OUT  CLUTTER  HAP  THRESHOLD  ATTENUATION 

C 

PRINT  927 

927  FORMAT <///»lX»7X» SFREQS r 1 OX ' ST ATTNS ' 1 OX ' STATTNDBS » / ) 

DO  14S  I-1'NPTS 

PRINT  928 * I * FREQ ( I ) * T ATTN  < I > » TATTNDB ( X > 

928  F0RMAT<1X'I2'3X'F7.2'E15.8'2X'E15.8) 

149  CONTINUE 

C 

C CHECK  IF  GROUND  CLUTTER  PLU8  N0I8E  ATTENUATION  VS. FREQUENCY  IS  TO  BE  PLOTTED 

- C 

199  IF < .NOT » IATTN)GO  TO  178 

C 

_ C COMPUTE  GROUND  CLUTTER  ATTENUATION  AND  CONVERT  TO  DB 

C 

DO  172  K-1'NPTS 

TEST-SIGCSGK  1 >SMUCL+MUNS8XGN8Q( 1 > 

IF < TEST . EQ . 0 . 0 ) TEST-1 . E-SO 

ATTN < K ) - < MUCLS8XGC8Q < K ) +HUNS8 XGNSQ  < K > > /TEST 

IF< ATTN(K) . EQ.O. ) ATTN<K)«1 .E-90 

ATTNDB(K>-10.0SAL0G10<ATTN(K>) 

IF<ATTNDB(K> .LT. <-100. 0) >ATTNDB<K>— 100.0 
172  CONTINUE 

_ C 

C COMPUTE  FREQUENCY  ARRAY 

C 

FACT-1. 0/<TSNPTS> 

DO  174  K-1.MPT8 
FREQ<K)«<K-1 )SFACT 
174  CONTINUE 

- C 

C PLOT  GROUND  CLUTTER  PLUS  N0I8E  ATTENUATION  VS. FREQUENCY 

C 

PRINT  929 

929  FORMAT <//* IX 'SQROUND  CLUTTER  ATTENUATION  V8.FREQ.PL0TS'/' 

1 1X.S8ET  UP  PAPER  FOR  0 TO  -80  DB.  AND  TYPE  QOS) 

READ  804 'XGO 

CALL  PLOT < FREQ 'ATTNDB'NPTS) 
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C PRINT  OUT  GROUND  CLUTTER  ATTENUATION 
C 

PRINT  930 

930  F0RMAT<///*lX*7X*«FREQ*»10X*t  ATTN*»10X*t  ATTNDB**/) 
DO  177  X-1*NPTS 

PRINT  928* X *FREQ( X > *ATTN(X > * ATTNDB < X > 

177  CONTINUE 

178  1F< .NOT. IFRESPIGO  TO  170 
DO  171  K-i*NPT8 
TEST«8I0N8Q( 1 > 

IF(TE8T.EG.0.0)TE8T«1 . E-50 
FRE8P ( K > -SXGNSO  < K > /TEST 
IF(FRE8P(K> .EQ.O. )FRE8P<K)-l.E-50 
FRE8PDB ( K > - 1 0 . t ALOQ 1 0 ( FRE8P ( K ) ) 
IF<FRE8PDB<K).LT.-100.)FRE8PD»(K)— 100. 

171  CONTINUE 

FACT-1 ./<T*NPT8> 

DO  170  K-l *NPT8 
FREQ(K)-<K-1)«FACT 
170  CONTINUE 
C 

C PLOT  OF  FILTER  RESPONSE 

C 

PRINT  183 

183  FORMAT <// * IX  * ^FREQUENCY  RESPONSE  08  TOT  UEL»> 

READ  804*180 

CALL  PLOT  < FREQ  » FRE8PDB  » NPT8  > 

C 

C CHECK  TO  SEE  XF  CFAR  TO  BE  CONTINUED 

C 

IF ( .NOT.M0RCFAR)STOP11 
C 

C COMPUTE  8IGNAL-T0-INTERFERENCE  08. TARGET  DOPPLER 

C 

170  DO  182  K-1*NPT8 

T I ( K ) -HUCL*6IQC8Q< K > + MUM0«8XQM8Q< K ) +MUN*8I8N8Q< K ) 

182  CONTINUE 
C*«***tt*t* 

C OUTPUT  TI(K) 

C********** 

C PRINT  931 

931  F0RMAT<//*1X*«  TI  ARRAY**/) 

C PRINT  902* (TKK> *K«1*NPT8> 

902  FORMAT (/* IX *0E13»0*/*10( IX *0E13»0*/>  > 

C 

CALL  T8CALE ( WGT8 * MTX  * MTXMDL  * NPT8 * T * 8CALET * XWOT  > 

C WRITE  0UT5 

C 

DO  185  K-1*NPT8 
TE8T-TKK) 

IF < TE8T . EQ . 0 . ) TE8T- 1 . E-50 
8TI < K > • < MUT*8CALET < K > > /TEST 
IF<8TX<K).EQ.0.>STX<K>-l*E-50 
8TIDB(K)-10.0*AL0G10(STX<K> > 

IF(STIDB<K> .OT • 100. 0)8TIDB(K >-100.0 
IF<8TIDB(K> .LT. < -100.0) )8TIDB(K >— 100. 0 
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183  CONTINUE 

XF<.N0T.8TIPL)G0  TO  4010 
F ACT-1. /<T*NPT8> 

DO  4020  K-1»NPTS 
FREQ (K)« (K-l )*F ACT 
4020  CONTINUE 

PRINT  1000 

1000  FORMAT!//* lXf*BXGNAL-TO-XNTERFERENCE  PLOT*) 

READ  804*100 

CALL  PLOT < FREQ *8TXDB*NPT8> 

4010  CONTINUE 

TEMPORARY  STOP 

180  CONTINUE 

IF(THMDL.EG»1)G0  TO  3000 
8T0P13 

3000  IF(RH00.EQ.1)G0  TO  5010 
ST0P15 

C CLUTTER  NAP  THRESHOLD  CONSTANT 
3010  XKH-10.**(-6./NCLHAP> 

XKH-NCLHAP/XKH-NCLHAP 
WRITE <6 *5020 )XKM 
5020  FORMAT ( //6X*CLUTTERNAP  FALSE  ALARM  C0N8TANT-*E20.8> 

C CFAR  CONSTANT  DETERMINATION 
XNRD2-NRCELLS/2. 

XKCFAR- 10 . **  < 6 . /XNRD2 > 

XKCFAR-XNRD2* ( XKCFAR- 1 . > 

WR I TE  ( 6 . 3050  ) XKCF  AR 

SOSO  FORMAT <6X*CFAR  FALSE  ALARM  C0NSTANT-*E20.8> 

C AOG  CLUTTERHAP  THRESHOLDS 
CALL  WBULL1(MUCL»A*8XG) 

C WRITE  0UT1 

DO  3061  K-l *NPT8 
PFASUM(K)-0. 

SUMPD<K)-0» 

5061  CONTINUE 

DO  3030  NW-1'NWW 
DO  3040  K-ltNPTS 

T I 1 < K ) -SIG ( NW ) *8IGCSQ ( K > +MUWO*SIOWSQ ( K > +MUN48XGN8Q < K ) 
TE8T-TIKK) 

IF (TEST .EQ.O. > TEST-1 ,E-30 
ST 1 1 < K > -MUTB8CALET ( K ) /TE8T 
CMT(K)-TI1<1)STATTN(K> 

3040  CONTINUE 
C AVG  CFAR  THRE8H0LD8 
DO  3060  K-1»NPT8 

CFAR1 (K)-(MUW1-MUW0)*8IGW8Q(K)+TI1<K) 

CFAR1 < K > -CFAR1 ( K > BXKCFAR 
CF AR2 < K ) -T I 1 < K ) + < HUW2-MUW0 > *8I0W8Q < K > 
CFAR2(K>-CFAR2<K> BXKCFAR 
3060  CONTINUE 

DO  3080  K-l » NPTS 
TE8T-TIKK) 

IF < TE8T . EQ . 0 . ) TEST- 1 . E-30 
BETA<K)-TATTN<K)*TI1 ( 1>/TE8T 
PFA(K)-<NCLNAP/(NCLMAP4XKM*BETA<K> > >**NCLHAP 
PF  A8UM  < K ) -PF  A < K > +PFA8UM  < K ) 
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5080  CONTINUE 
6001  DO  6000  K-1*NPTB 

TTE8T < K > -AMAX 1 ( CNT < K > » CFAR1 < K > * CF AR2 < K ) > 

IF<K.E0.1)G0  TO  6010 
IF (TTE8T<K) .EO.CHT(K) )00  TO  6010 
IF(TTE8T (K> .EQ.CFARKK)  >00  TO  6030 
NTHOLD-3 
00  TO  6100 
6010  NTHOLD-1 

00  TO  6100 
6030  NTHOLD-2 

6100  GO  T0<6110*6120*6130) *NTHOLD 
6110  XKP-XKM 

NP-NCLMAP 
E<K)-BETA<K) 

00  TO  6140 
6120  NP-XNRD2 

XKP-XKCFAR 
TE8T-TIKK) 

IF  < TEST . EO . 0 . ) TEST-1 . E-30 
E < K ) -CF AR 1 < K > / < XKCF AR*TE8T ) 

00  TO  6140 
6130  XKP-XKCFAR 
NP-XNRD2 
TE8T-TIKK) 

IF < TEST . EO . 0 .) TE8T-1 . E-50 
E < K > -CF AR2 < K ) / ( XKP8TE8T > 

6140  PD<K)-1 •/< 1 » + <XKP*E<K>  >/<NP*< 1 •♦8TI1<K) > ) )8*NP 
8UMPD ( K ) -8UHPD  C K ) +PD  < K ) 

6000  CONTINUE 

C IF<NU»LE»2)WRITE  0UT2 
5030  CONTINUE 
5090  DO  6160  K-lfNPT8 

FRE0<K)-(K-1.)*1./<T*NPT8) 

6160  CONTINUE 

DO  6150  K-l rNPTS 
PFA8UM ( K > -PFA8UH < K > /NWW 
SUHPD < K ) -SUNPD ( K > /NUN 
6150  CONTINUE 

IF < .NOT *PRFA)00  TO  6200 
C PROBABILITY  OF  FAL8E  ALARM  V8  FREQUENCY 
WRITE(6f 6169> 

6169  FORMAT (6X8PR0B.  OF  FALSE  ALARM  US  FREQUENCY*) 

READ  804 f 100 

CALL  PLOT ( FREQ »PFASUM»NPT8> 

WRITE<6r 6170) 

6170  FORMAT < //6X*FREQUENCY*10X*PR0B  OF  FALSE  ALARM*//) 

MRITE(6»6310) (FREQ(K) rPFASUM(K) » K-l rNPTS) 

C PROBABILITY  OF  DETECTION 
6200  CONTINUE 

IF ( «NOT »PRED)GO  TO  7001 
MRITE<6r6210) 

6210  FORMAT (6X*0RAPH  OF  DETECTION  PROBABILITY*) 

READ  804*100 

CALL  PLOT ( FREQ *8UMPD*NPT8) 

URITEf  6*6300) 

6300  FORMAT<//6X*FREGlUENCY*10X*DETECTION  PROBABILITY*//) 
WRITE<6»6310) (FREQ(K) »8UMPD(K) *K-1 *NPTS) 
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6310  FORMAT < 1X*E20.8*4X*E20»8) 

7001  IF ( »H0T » IDEALCF)G0  TO  7000 
YB-13. 81551056 
DO  7090  K-1*NPTS 
SUMPDOO-O. 

7090  CONTINUE 

DO  8000  NW-1*NWW 
DO  8010  K-1*NPT8 

T 1 1 < K )-8I0<  MU > *8 IQC8Q ( K > +MUW0B8 1 GWSQ ( K > +MUN*8IGN8Q< K ) 
TEST-TIKK) 

IF (TEST «EQ«0« )TE8T"1 *E~50 
STI(K) -MUTB8C ALET <K>/TE8T 
PD(K)-EXPC-YB/C1.+STI(K>>) 

80HPD ( K ) -8UMPD < K ) + PD < K > 

8010  CONTINUE 
8000  CONTINUE 

DO  8020  K-l f NPT8 
SUMPD ( K ) -8UHPD ( K > /NWW 
8020  CONTINUE 

WRITE (6 *8030) 

8030  FORMAT (6XBPR0BABILITY  OF  DETECT ION t IDEAL  CFARB) 

READ  804*180 

CALL  PLOT < FREQ *8UMPD*NPTS) 

WRITE (6*6300) 

WRITE(6*6310) (FREQ(K) *8UMPD(K) *K-1 *NPT8) 

7000  CONTINUE 
CALL  EXIT 
END 

C*****««t*******««*******«***«*»«*****t«t*»*M«****«M«t*«*t****** 
C*«*************«******* 8t8*********8**8*****8****88t*««B*B**** 
SUBROUTINE  T6CALE ( WOTS  * HTI * HTINDL *NPT8 * T * 8CALET * I WOT ) 


LOGICAL  IWQTrhTI 

COMPLEX  X * XTEMP  * SUM » AVQX » XBAR 

DIMENSION  WOTS ( 64 ) * SCALET ( 64 ) * X ( 66 ) » XTEMP (64) *XBAR(64) r FREQ (64) 


PI-3.14159265 


C COMPUTE  FREQUENCY  ARRAY 
C 

FACT-1 »0/(TBNPT8) 

DO  10  K-1*NPTS 
FREQ (K)-(K-l) BFACT 
10  CONTINUE 

C 

C COMPUTE  X(L>  »L-1 *NPT8P2  FOR  EACH  K 
C 

NPT8P2-NPT8+2 

C 

C BEGINNING  OF  K LOOP 
C 
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DO  15  K-lfNPTS 
C 

DO  20  L-1,NPT8P2 
AR0-2.0*PI*FREQ(K)*<L-1)*T 
X<L)-CMPLX<COB< ARQ) »8IN< ARG> ) 

20  CONTINUE 
C*****«**«t 

C OUTPUT  X<L>  FOR  EACH  K 
C***««*t*»t 
C PRINT  900?K 

900  FORMAT <//r*  X(L>  FOR  K-*'I2'/> 

C PRINT  902» <X<L).L-1»NPT8P2> 

902  F0RMAT</tlX'2E13*6f2Xr2E13.6'/f31<lX»2E13.6r2X'2Ei3.6'/>> 
C 

C CHECK  IF  MTI  TO  BE  DONE 
C 

IF ( .NOT. MTI >00  TO  30 
C 

C IF  MTI  TO  BE  DONEr  MODEL  *1  OR  MODEL  92  T 
C 

IF(MTIMDL.EQ.2)00  TO  25 
C 

C MTI  MODEL  «1 
C 

DO  22  L-l »NPT8 

XTEMP<L)-X(L>-2.0*X<L+1 )+X(L+2) 

22  CONTINUE 

DO  23  L-1.NPT8 
X<L)-XTEMP<L> 

23  CONTINUE 
C«****«*«*« 

C OUTPUT  X(L)  FOR  EACH  K AFTER  MTI  MODEL  91 
C***»**t*tt 
C PRINT  903*K 

903  F0RMAT<//»1X>*  X<L>  AFTER  MTI  MODEL  91  FOR  K»*tI2»/> 

C PRINT  902»  <X<L) »L-lrNPTS> 

C 

00  TO  30 
C 

C MTI  MODEL  92 
C 

25  SUM-0.0 

DO  35  L-1»NPT8 
SUM-8UM+X ( L > 

35  CONTINUE 

AUQX-8UM/NPT8 
DO  37  L-1»NPT8 
X<L)-X<L)-AVGX 
37  CONTINUE 
C***««*t**» 

C OUTPUT  X(L)  FOR  EACH  K AFTER  MTI  MODEL  92 
C*«*»t*«*tt 
C PRINT  904»K 

904  FORMATf // * IX » * X<L>  AFTER  MTI  MODEL  92  FOR  K«*»I2»/> 

C PRINT  902'<X<L)'L-1'NPTS> 

C 
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C APPLY  WEIGHTS  IF  REQUIRED 
C 

30  IF  < *N0T « IWGT)G0  TO  40 
DO  39  L-1pNPTS 
X<L)-X<L>*WGT8<L> 

39  CONTINUE 

C»****«***« 

C OUTPUT  X(L)  FOR  EACH  K AFTER  WEIGHTING 

C********** 


PRINT  90S pK 
905  F0RMAT(//p1Xp* 


X<L)  AFTER  WEIGHTING  FOR  K-*pI2p/> 


PRINT  902f  <X(L) pL«1pNPTS) 


C COMPUTE  XBAR(K) 

C 

40  SUM-0.0 

DO  45  L-1pNPTS 

ARG» < 2 . 0*P I > /NPTS* < < L- 1 ) * < K- 1 ) ) 

SUM-SUM+X ( L ) SCMPLX ( COS ( ARG  > p -S IN < ARG  > > 

45  CONTINUE 
XBAROO-SUM 
C********** 

C OUTPUT  XBAR(K) 

C*«***t**«* 

C PRINT  906»Kf XBAR(K) 

906  F0RMAT(//»lX»*K-*»I2»2Xr*XBAR(K)-*f2E13.6f/) 

C 

C DEFINE  OUTGOING  SCALE T ARRAY 

C 

SCALET ( K ) - ( REAL < XBAR <K>  > > **24 ( AIHAG ( XBAR < K ) > > **2 
C********** 

C OUTPUT  SCALET (K) 

C********** 

C PRINT  907rK»8CALET(K> 

907  F0RNAT(//»lX»*K-*»I2r2Xr*SCALET(K>-*rE13.6r/) 

C 

C END  OF  K LOOP 
C 

15  CONTINUE 
C 

RETURN 

END 

SUBROUTINE  PLOT(X»YtN> 

LOGICAL  IPROB 
DIMENSION  X<64>  pY<64> 

DIMENSION  TEMP (64) 

NAMELIST/INP1/YMIN»YMAX 
NAMELI8T/INP2/IPR0BfXMIN»XMAX 
NAMEL I ST/INP3/XMIN  p XMAX » YMIN  » YMAX 
DO  1 I-1pN 
TEMP< I >«Y< I > 

1 CONTINUE 
XMIN-X(l) 

XMAX-X(l) 

YMIN-TEMP(l) 

YMAX»TEMP< 1 ) 

DO  10  I-1pN 
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DIFXl»X< I >-XMIN 

DIFX2-X<  X )-XMAX  i 

DIFY1-TEMP<I)-YMIN  | 

DIFY2-TEMP<X>-YMAX  1 

IF<DIFX1.LT.0.)XMIN-X<I>  j 

IF (DIFX2.GT . 0. )XMAX-X< I ) 

IF  < DIFY1 .LT.O. ) YMIN-TEMP  < I > 

IF < DIF Y2 . BT . 0 . ) YHAX-TEHP < I > 

10  CONTINUE 
WRrTe  I MP3 
WRITE(6*11 ) 

11  FORMAT <6X*INPUT  INP2t IPROB*XMIN*XMAX«> 

READ  INP2 

IF  < .NOT • XPR0B)Q0  TO  12 
CALL  P8CALE<TEMP*N*YMXN*YMAX> 

00  TO  13 
PRINT  INP1 

12  WRXTE<6*20> 

20  F0RMAT<6X*INP  NAMELIST  INP1 I YMIN»YMAX*> 

READ  XNP1 

13  PRINT  924 

924  FORMAT < IX **PLTL*> 

NN-N/2.+1. 

DO  30  1-1 »NN 

I X- < X ( X ) -XM XN > *9999 . 0/ ( XMAX-XHI N ) 

I Y- < TEMP < I ) - YMXN ) *9999 . 0/ < YMAX- YM IN > 

PRINT  92S»IX*XY 

925  FORMAT (2< IX* 14)  > 

30  CONTINUE 

PRINT  926 

926  FORMAT < IX »*PLTT*> 

RETURN 

END 

SUBROUTINE  WBULL1 (MUCL*A*8IG) 

REAL  MUCL 

DIMENSION  SXGUOOO) 

DIMENSION  R<1000> 

DO  10  1-1*1000 
R(I)-RANF(0> 

SIQ( I )*( ALOO( 1»/<1.-R(I)>> >««A 
SXO<  X >-MUCL*8I0( I >/2. 

10  CONTINUE 
SUM-0. 

DO  20  X-l»1000 
SUM-SUM+8IG< I > 

20  CONTINUE 

8UM-SUM/1000. 

DO  30  1-1*1000 
810 < I ) - MUCL *810 ( I > /SUM 
30  CONTINUE 

WRITE <6*50) SUM 
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SO  FORMAT <6X*MEAN  OF  UN8CALED  CLUTTER  CROSS  SECTION  SEQ  I8*E20.8) 
RETURN 
END 

SUBROUT I NE  PSC ALE < Y » NPT8 f YM IN , YMAX > 

DIMENSION  Y<64> 

F I < P > -<  A0+A1*P+A2« < P**2 > ) / < 1 . +B1BP+B2* < P**2 > +B3* ( P**3 > > 

AO-2.313317 

A 1-. 802853 

A2-.010328 

Bl- 1.432788 

B2-.1892&9 

B3-.001308 

AS-. 05 

PS-SORT ( ALOO  < 1 ♦ / < A8««2  > ) ) 

YMIN-FI (PS) 

YMIN-PS-YMIN 
YMIN— YMIN 
AL-.0001 

PL-SQRT  < ALOO  < 1 . / < AL**2 ) > ) 

YMAX-FKPL) 

YMAX-PL-YMAX 

DO  10  I-l »NPTS 

I F ( Y ( I > . LE . 0 • 5 ) GO  TO  20 

AA-1.-1.E-10 

IF<Y(I).EQ,1.)Y(I )-AA 

P-SORT  < ALOG  < 1 . / < ( 1 . -Y ( I ) >**2 ) > ) 

Y(I)-FKP) 

Y< I )-P-Y< I > 

00  TO  10 

20  IF< Y< I ) . EQ.O. ) Y< I >— 1 .E-10 
P-SQRT ( ALOO < 1 , / < Y < I ) **2 ) ) ) 

Y < I >-FI <P> 

Y(I)-P-Y<I> 

Y < I )— Y < I ) 

10  CONTINUE 

RETURN 

END 
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j - APPENDIX  B 

USE  OF  EXISTING  FAST  FOURIER  TRANSFORM  DIGITAL 
! “ PROCESSOR  (FFTDP)  FOR  EAR 


A block  diagram  representation  of  the  FFTDP  is  depicted  below.  As 
shown  it  consists  of  a pair  of  9-bit  A/D  converters,  a time  weighting 
circuit,  an  input  storage  buffer,  and  a complex  adder/subtractor/multiplier 
arranged  in  an  FFT  butterfly  configuration  which  outputs  to  a post  processor. 


Parameters  of  the  FFTDP  are: 

1.  9-bit  A/D  with  a 2.4  usee  (integrate  and  dump)  sample  rate. 

2.  Input  buffer  of  512  words  of  48  bits. 

3.  (Cosine)^  time  weighting  circuits. 

4.  64  point  FFT  (complex  adder/complex  subtractor/complex 
multiplier  butterfly  configuration). 

5.  30  range  gates  plus  1 test  gate. 

6.  Fixed  4.608  ms  time  to  process  a 64  point  batch. 

The  organization  of  the  input  buffer  is  depicted  below;  it  is  made  up 
of  512  words  of  48  bits.  Each  48-bit  word  is  made  up  of  4 data  words. 

Hence,  the  FFT  arithmetic  unit  carries  only  12  bits  with  a scaling  capability 
when  the  arithmetic  overflows  the  12  bits. 
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Address 


S!R1 

S1R2 

S33R1 

S33R2 

S1R3 

S1R4 

S33R3 

S33R4 

S1R5 

S33R5 

S1R29 

S1R30 

S33R29 

S33R30 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

S32R29 

S32R30 

S64R29 

S64R30 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

Spare 

R = Range  Gate 

Memory  Organization 

The  relationships  within  the  buffer/sample/FFT  process  are  shown 
pictorially  below: 


4.  608  ms 


m 72  US  ►j 

2.  4 us/RG  x 30  s 72  us 


mrs  = 13'888-8 


process  time  = 4.608  ms 

time  of  1 Butterfly  = Ta"””””""  = = 800  ns 

Tlog264 

To  use  this  processor  with  the  EAR,  a new  synchronizer  must  be 
designed  as  well  as  another  input  buffer  suitable  for  reformatting  the  input 
data  into  blocks  of  30  range  gates.  The  64  point  input  data  is  collected  by 
the  EAR  in  64  x 200  = 12,800  usee  at  its  nominal  5 kHz  PRF.  The  FFTDP 
processes  a batch  of  30  range  gates  worth  of  data  in  4608  msec.  Therefore 
12,800/4,608  = 2.777  batches  can  be  processed.  This  must  be  rounded 
downward  to  the  nearest  integer,  so  that  2 x 30  ■ 60  range  gates  can  be 
handled  by  the  processor  at  the  5 kHz  PRF.  It  is  probably  sensible  to 
lower  the  maximum  PRF  to  4629  Hz.  The  data  collection  time  then  goes 
up  to  64  x 216  = 13,826  msec. , which  will  then  handle  exactly  three  batches, 
for  a total  of  90  range  gates. 

With  such  a small  number  of  range  gates,  the  kind  of  CFAR 
processing  to  be  used  is  open  to  question.  A sliding-window  cell-averaging 
CFAR  as  discussed  in  the  text  may  not  make  too  much  sense,  since  end 
effects  will  tend  to  dominate. 


> 


f 

i 


! 


B 


# 


DISTRIBUTION 


Oefense  Docunentatlon  Center 
Cameron  Station 
Alexandria,  Virginia  22314 

Commander 

U.S.  Army  Materiel  Development  and  Readiness  Conmand 
ATTN:  DRCCG 
DRCRD 

5001  Elsenhower  Avenue 
Alexandria,  Virginia  22304 

Commander 

Ballistic  Missile  Defense  Systems  Command 
ATTN:  BMDSC-HR,  Mr.  Stevenson 
P.0.  Box  1500 

Huntsville,  Alabama  35807 
Director 

Ballistic  Missile  Defense  Advanced  Technology  Center 
ATTN:  ATC-R,  Mr.  Carlson 
P.0.  Box  1500 
Huntsville,  Alabama  35807 

Commander 

US  Army  Electronics  Command 
ATTN:  DRSEL,  Mr.  Flshblen 
DRCPM-MALR 

Fort  Monmouth,  New  Jersey 

DRCPM-MDE,  Mr.  Wood 
DRCPM-HAE , Mr.  Ams 
DRCPM-CFE,  Mr.  David 
DRCPM-SHO , Mr.  Bishop 
DRSMI-R,  Dr.  McDaniel 
-R,  Mr.  Fagan 
-R,  Dr.  Kobler 
-RE,  Mr.  Llndberg^ 

-RE,  Mr.  Pittman 
-REO,  Mr.  Currie 
-RER,  Mr.  Cash 
-REG,  Dr.  Burl  age 
-RE,  Record  Set 
-R6,  Mr.  Huff 
-RD,  Dr.  McCorkle 
-RBL 

-RP,  Dr.  Jackson 


1 

No.  of  Copies 
12 


1 

1 


1 


2 


50 


5 

1 


